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Abstract 

We study a new class of rank two sub-Riemannian manifolds encompassing Rieman- 
nian manifolds, CR manifolds with vanishing Webster- Tanaka torsion, orthonormal bundles 
over Riemannian manifolds, and graded nilpotent Lie groups of step two. These manifolds 
admit a canonical horizontal connection and a canonical sub-Laplacian. We construct on 
these manifolds an analogue of the Riemannian Ricci tensor and prove Bochner type formu- 
las for the sub-Laplacian. As a consequence, it is possible to formulate on these spaces a 
sub-Riemannian analogue of the so-called curvature dimension inequality. Sub-Riemannian 
manifolds for which this inequality is satisfied are shown to share many properties in common 
with Riemannian manifolds whose Ricci curvature is bounded from below. 
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1 Introduction 

A sub-Riemannian manifold is a smooth Riemannian manifold M equipped with a fiber inner 
product •) on the tangent bundle TM and a non-holonomic, or bracket generating, sub- 
bundle 7i C TM. This means that if we denote by L(J-L) the Lie algebra of the vector fields 
generated by the global C°° sections of TL, then span{X(x) | X G L{TL)} = Tx{M) for every 
X G M. A piecewise smooth curve 7 : [a, 5] ^ M is called admissible, or horizontal, if it is 
tangent to Ti, i.e. if 7'(t) G whenever 7'(t) is defined. The horizontal length of 7 is 

defined as in Riemannian geometry 

^h(7)= f gR[i{t),i{t))dt. 

J a 

Denoting by y) the collection of all horizontal curves joining x,y € M, one defines a distance 
d{x,y) between x and y by minimizing on the length of all 7 G TC{x,y), i.e. 

d{x,y) = inf £n{l)- 

Such distance was introduced by Caratheodory in his seminal paper [lS| on formalization of the 
classical thermodynamics. In such framework horizontal curves correspond loosely speaking to 
adiabatic processes. 

To be precise, in [16] the question of whether d{x, y) be a true distance was left open. This 
question was answered by the fundamental connectivity theorem of Chow [20j and Rashevsky ^54j 
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which states that if M is connected and 7i is bracket generating, then 7i{x, y) ^ for every x,y £ 
M. As a consequence, d{x, y) is finite and therefore it is a true distance. Such metric is nowadays 
known as the control, or Carnot-Caratheodory distance on M (after Gromov, Lafontaine and 
Pansu [Ml, see also [Si], [SH]). Besides the cited references the reader should consult E. Cartan's 
pioneering address [TB] at the Bologna International Congress of Mathematicians in 1928, as well 
as the articles by Rayner [55] (where sub-Riemannian manifolds are called parabolic spaces), 
by Mitchell [50] and Strichartz [6l], see also [65]. One should also consult the monographs [lOj . 
[SI], [I] and [8]. 

We note that when TC = TM, then the distance d{x, y) is simply the Riemannian distance associ- 
ated with the inner product 5_r(-, •), and thus sub-Riemannian manifolds encompass Riemannian 
ones. However, some aspects of the geometry of sub-Riemannian manifolds are considerably less 
regular than their Riemannian ancestors. Some of the major differences between the two ge- 
ometries are the following: 

1. The Hausdorff dimension of the metric space (M, d) is usually greater than the manifold 
dimension; 

2. The exponential map defined by the geodesies of the metric space (M, d) is in general not 
a local diffeomorphism in a neighborhood of the point at which it is based (see |"55j ) ; 

3. The space of horizontal paths joining two fixed points may have singularities (the so-called 
abnormal geodesies, see [ST]). 

During the last two decades there have been several advances in the study of sub-Riemannian 
spaces and the closely connected theory of sub-elliptic pde's, see |62], [28], [58], [26], [52j, [61] . 
[39] . [38] . [27], [H], [12], [7l]. However, these developments are of a local nature. As a conse- 
quence, the theory presently lacks a body of results which, similarly to the Riemannian case, 
connect properties of solutions of the relevant pde's to the geometry of the ambient manifold. 

The purpose of this paper is to begin a program aimed at filling this gap. Precisely, in a sub- 
Riemannian manifold of rank two we introduce a new notion of Ricci curvature tensor, and with 
such notion we obtain various results which parallel those of Riemannian manifolds with Ricci 
tensor bounded from below. 

To put our work in the proper perspective we recall that in Riemannian geometry a first point 
of view on the Ricci tensor is to understand it as a measure of volume distortion in geodesic 
normal coordinates. This point of view has recently led several authors ( Lott-Villani [49j, 
Sturm [66], [67], Ollivier [53]) to introduce on metric spaces more general than the Riemannian 
ones a suitable notion of Ricci curvature based on the theory of optimal transport (see [72]). 
However, as pointed out in [40], this notion of Ricci curvature bound may not be applied in a sub- 
Riemannian framework. In sub-Riemannian geometry, the quantification of volume distortion 
properties is particularly difficult to handle because of the singular nature of the exponential 
map. To the authors' best knowledge, the only significant results in this direction are at the 
moment known only in the special case of three-dimensional contact manifolds (see [59], [IQ] 
and [2]). 
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There is a second point of view on the Riemannian Ricci tensor and this is the one that we have 
adopted in the present paper. On a Riemannian manifold, there is a canonical second order 
differential operator: the Laplace-Beltrami operator A. Properties of this operator and of the 
associated heat flow are intimately related to the geometry and the topology of the underlying 
manifold. One of the cornerstones of the interplay between the analysis of the Laplace-Beltrami 
operator and geometry is given by the celebrated Bochner formula: 

A(|V/|2) =2||vV||' + 2< V/,V(A/) >+2 Ric(V/,V/), 

This is where the Ricci curvature tensor appears in the study of the Laplace-Beltrami operator, 
and it is then seen that a lower bound assumption on the Ricci curvature is equivalent to a 
coercivity property of a canonical bilinear differential form associated to A. More precisely, 
associated with A are the two following differential bilinear forms on smooth functions /, g : 

M ^ R, 

r{f,g) = ^(A(/g) - fAg - gAf) = {Vf,Vg), 

and 

T2if, g) = ^ [AT{f, g) - T{f, Ag) - T{g, A/)] . 
As an application of the Bochner's formula, which we can re-write as 

Ar(/, /) = 2| I VVl I' + 2r(/, A/) + 2 Ric(V/, V/), 

one obtains 

r2(/,/) = ||v2/i + Ric(v/,v/). 

With the aid of Schwarz inequality, which gives 

llVVlli > the assumption that the 

Riemannian Ricci tensor on M is bounded from below by /9 G R, translates then into the so- 
called curvature-dimension inequality: 

r2(/,/)>i(A/)2 + pr(/,/). (Li) 

The inequality (jl.ip perfectly captures the notion of Ricci curvature lower bound, see example 
15.21 below. In the hands of D. Bakry, M. Ledoux and their co-authors it has proven to be a 
powerful tool in recovering most of the well-known theorems which, in Riemannian geometry, 
are obtained under the assumption that the Ricci curvature be bounded from below (see for 
instance [3], [45], [47]). 

The purpose of this work is to generalize this point of view in a sub-Riemannian setting. More 
precisely, we assume that on M we are given smooth vector fields Xi,...,X(i which generate 
the horizontal subbundle 7i. The commutators [Xi^Xj] are supposed to satisfy some structural 
assumptions, see (|2.3p and (j2.4|) below, with Z^n^m^n = being the non-horizontal, or 

vertical directions. Setting 

^.{x) = span{Xi(x), V{x) = s\ian{Zmn{x) | 1 < m < re < [)}, 
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we assume that 

T^M = n{x)®V{x), xGM. 
We also assume that Ti. is bracket-generating of rank two, i.e. 

TM = span{X,(x), [Xj,Xk]ix)}, xGM. 

The first step wiU be to equip M with a canonical subelliptic operator 

d 
i=l 

which shall play the role of the Laplace-Beltrami operator. We then equip M with a degenerate 
metric tensor g{-, ■) for which Xd{x)} is orthonormal at every x G M, the spaces 

and V(x) are orthogonal, and such that ^Iv = 0. In such a manifold we introduce a canonical 
connection V generalizing the Levi-Civita connection. One of the fundamental assumptions in 
this work is that the torsion be vertical. Besides the Riemannian case, basic examples covered by 
our framework are CR manifolds with vanishing Webster- Tanaka torsion (Sasakian manifolds), 
orthonormal bundles over Riemannian manifolds, and graded nilpotent Lie groups of step two. 

An essential tool in our analysis are two sub- Riemannian Bochner formulas for L, one in the 
horizontal direction and one in the vertical one. Such formulas are established in section HI We 
mention here that their technical complexity is the main raison d'etre of the above rank two 
assumption (in this regard, one should see the closing comments at the end of this introduction). 
By means of these formulas we identify a tensor TZ which plays the role of the Riemannian Ricci 
tensor. We show in Proposition 13.11 that for every / G C°°(M) one has 

d 

nfJ) = E ^^<Xe,Xk)XefXkf - {{V x,T){X,,Xk)f){Xkf) + - {T{X,,Xk)ff . 

l,k=l 

where T is the torsion of the canonical connection V and Ric its Ricci curvature. In fact, in the 
Riemannian case, one has 

7^(/,/) = Ric(V/,V/), /gC°°(M). 

A lower bound assumption on the tensor IZ will translate into a generalized curvature dimension 
inequality for L that writes in the form: 

r2(/, /) + i^rf (/, /) > i(L/)2 + - ^) r(/, /) + p^v^U. f), ^>o, (i-2) 

where and are two bilinear differential forms on the vertical subbundle. The parameters 
of this inequality are pi,p2,K, and a recurrent theme of our work is that all estimates solely 
depend in an explicit quantitative way on these parameters. As shown in the following summary 
of the main results, the parameter pi is of particular importance: 

1. Dodziuk-Yau type theorem: If the inequality ()1.2p holds for some constants pi £ 
M, P2 > 0, K > 0, then the heat semigroup Pt = e*-^ is stochastically complete and bounded 
solutions of the heat equation are characterized by their initial condition; 
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2. Li-Yau type inequality: If the inequality (|1.2p holds for some constants pi £ M, p2 > 
0, K > 0, then the heat semigroup Pt satisfies a Li-Yau type gradient estimate. Exploit- 
ing the latter, we prove that the heat kernel p{x, y, t) of Pt satisfies a uniform Harnack 
inequality and a Gaussian upper bound estimate; 

3. Yau-Liouville type theorem: If the inequality (jl.2p holds for some constants pi > 
0, p2 > 0, K > 0, then there is no non constant bounded harmonic function; 

4. Volume and isoperimetry: If the inequality (jl.2p holds for some constants pi > 0, P2 > 
0, K > 0, then the volume growth of the geodesic balls is at most polynomial. Moreover, 
under the assumption that the volume of balls is at least polynomial, we obtain a global 
isoperimetric inequality; 

5. Myers type theorem: If the inequality (|1.2p holds for some constants pi > 0,p2 > 
0, K > 0, then the metric space (M, d) is compact in the metric topology with a Hausdorff 
dimension less than d (^1 + and we have 



diam M < 2\/3- 




6. Poincare inequality: If the inequality (jl.2p holds for some constants pi > 0, P2 > 
0, K > 0, then the following inequality is satisfied 

inf / \f-c\df,<6d(i + —),[^^[ vnr)dp. 
cmju^-' V 2p2j y dpiP2 Jm^ ^ 



7. Lichnerowicz type estimate: If the inequality (II. 2p holds for some constants pi > 
0, /92 > 0, K > 0, then the first non zero eigenvalue Ai of —L satisfies the estimate 



^P2 + K- 



The basic source of the above listed results is the study of monotone entropy type functionals 
of the heat semigroup. More precisely, generalizing [5], [7] and [9], one of the most important 
observations is that in our framework, the inequality (jl.2p implies that for the functionals 



^i{t) = Pt{{PT-tf)T{\nPT-tf)) , 

^2{t) = Pt {{Pr-t f (In Pr-tf)) ■ 
we have the following differential inequality 

- — >--^LPTf + -^PTf, 
2p2 / dp2 dp2 

where b is any smooth, positive and decreasing function on the time interval [0, T] and 

d /b" K b' 



b" Kb' \ 
V P2 J 
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Depending on the value of pi, a careful choice of the function b leads then to a generalized 
Li-Yau type inequality from which it is possible to deduce the above results. 

Some final comments are in order. We have mentioned above that the choice of working with 
sub-Riemannian manifolds of rank two is closely connected with the complexity of the Bochner 
type formulas in section [H On one hand, the rank two setting is rich enough to encompass at 
one time the case of Riemannian and CR (Sasakian) manifolds. On the other hand, following the 
program in section HI for manifolds of arbitrary rank r one would need to establish r Bochner type 
formulas (of decreasing complexity, the most difficult one being the horizontal one). Whereas 
it would be desirable to treat sub-Riemannian manifolds of arbitrary rank, we have felt that 
the increased technical difficulties connected with this endeavor would distract from the main 
ideas, and have consequently decided to defer the treatment of manifolds of rank >3 to a future 
study. We should also mention that most of the constants appearing in the main results in this 
paper are not optimal. This can be seen by considering the special setting of graded nilpotent 
Lie groups where one can use the underlying non-isotropic dilations to obtain sharper constants. 
Finally, we mention that our methods strongly rely on the assumption that the torsion of the 
canonical connection is vertical. 

In closing we mention that a pseudo-hermitian version of the Bonnet-Myers theorem for contact 
manifolds of dimension three was proved by Rumin, see Theorem 16 in [59]. We thank S. 
Webster for bringing this to our attention. 

2 The framework and its geometrical interpretation 
2.1 Prelimimaries 

Henceforth in this paper, M will be a smooth connected Riemannian manifold. We assume that 
Xi, ...,Xci are given smooth vector fields on M satisfying the following commutation relations: 

d f) 

[X„X,] = ^ 4.x, + ^ 7r^-„, (2.3) 

1=1 m,n=l 
d 

[Xi, Z-inn] = ^ df^n^e, (2.4) 
e=i 

for some smooth vector fields {Zmn}i<m,n<t^ and smooth functions ufj, 'y^^ and By 
convention Z^n = —Znm, ^fj = and 7^*^" = —7^" • We will assume that 

^imn = -^}mn^ i, £ = 1, d, and m, 71 = 1, f). (2.5) 

Note that (|2.5p implies = for i = 1, .., d, and m, n = 1, f). The assumption (|2.5p plays 
a pervasive role in the results of this paper. 

We shall moreover assume that the vector fields Xj's satisfy Hormander's finite rank condition 
of step two [37]: for every x G M, 

span{Xi{x), [Xj,Xk]{x), I < i < d,l < j < k < d} = T^M. 
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We denote by 

n{x) = span{Xi{x), ...,Xd{x)} , x G M, 
the set of horizontal directions at x, and by 

V(x) = span{Zmn{x), 1 < m < n < I)} , 

that of vertical directions. We shall assume that 

dim'H(x) = d, dimV(x) = ^ — ^, 

at each point x G M, and that 

n{x) © V(x) = T^M. 

The dimension of M is therefore d + ^ but such number will never explicitly appear in the 

results in this paper. We indicate with Ti = Ua'eM'^(^)' ~ ^xeM^^i^)^ respectively the 

horizontal and vertical subbundles of TM. 

Our goal is to study the second order subelliptic operator L given by 

d 



L = Y,Xf + Xo, (2.6) 



i=l 

where 

d 

Xo = - ^ u;1,X,. (2.7) 

i,k=l 

We will assume that with respect to the Riemannian measure of M, 

^mn — ~Zmn-, L* = L, (2.8) 

where denotes the formal adjoint of Zmn and L* the formal adjoint of L. We note explicitly 
that (p^ means that, if we set for every ^,■0 S Co°(M) 



< (j),ip >= / (pTpdfj., 
Jm 

then one has 

< Z^n4>, >= - < 0, Zmni^ > , < Lc/), ^ > = < (l),L'4)> . (2.9) 

We recall that, thanks to Hormander's hypoellipticity theorem [37j . distributional solutions to 
Lf = in M are C°° functions. We stress that, thanks to the second identity in the assumption 
()2.8p . the operator L can be realized as 



L = -Y,XtX,. (2.10) 



i=l 
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We also assume that M is endowed with a Levi-Civita connection with respect to which the 
Laplace-Beltrami operator is given by 

d 

A = ^Xf+ Yl ^L + ^o. (2.11) 

i=l l<m<n<i) 

Let us now give some examples that, at least locally, fit into the previous framework and that 
constitute a basic motivation for our study. 

Example 2.1 [Laplace-Beltrami operator on a Riemannian manifold] Let (M, g) be a d- dimensional 
connected Riemannian manifold with Levi-Civita connection V. Let Xi^...,Xd he a local or- 
thonormal frame around a point xq E M. In that case, we have 

d 

[x,,x,] = Vx^x, - Vx,x, = ^ (rf^. - r^^,) x^ 

k=l 

where T^j are the Christojjel symbols of the Levi- Civita connection. Thus, in this particular case 
we have 7[^" = in ()2.3p . The Laplace-Beltrami operator on M reads 

d 

A = J2xf + Xo, 

i=l 

where 

d 

^0 = - 'Y (j^ik - ^ki) ^i- 

i,k=l 

Finally, we can observe that L is symmetric with respect to the Riemannian measure on M. 

Example 2.2 [Graded nilpotent Lie groups of step two] Let G be a connected and simply con- 
nected nilpotent Lie group of step two. This means that its Lie algebra can be written as 
Q = Vi (B V2, where [Vi,Vi] = V2, and [Vi,V2] = {0}. Let Lx{y) = xy be the operator of 
left-translation on G, and indicate with dL^ its differential. If ei, ...^ed is an orthonormal basis 
of Vi, we indicate with Xi, ...^Xd, where Xj{x) = dLx{ej), the corresponding system of left- 
invariant vector fields on G. We assume that G is endowed with a left-invariant Riemannian 
inner product with respect to which {Xi^ ...,Xci} constitutes a global orthonormal frame. In this 
framework, we see that (j2.3p holds with 

Jij = 0, 

mn ^ _^mn ^1 rn ^ r> 
imn Inm 2 ' ni- -/- ii', 

lij^ = 0, otherwise, 
5^ =0. 

In view of (j2.10p we see that L = Yli=i-Xi ^-^ sub-Laplacian associated with Xi, X^, see 
f28^ . 164^ . In this case L is symmetric with respect to the bi-invariant Haar measure on M. 
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Example 2.3 [Horizontal Bochner Laplace operator] Let {M.,g) be a d- dimensional connected 
smooth Riemannian manifold endowed with the Levi-Civita connection. Let us consider the 
orthonormal frame bundle O (M) over M. For each x we can define a horizontal vector 

field Hx on O (M) by the property that at each point u G 0{M), Hx{u) is the horizontal lift of 
u{x) from u. // (ei, e^) is the canonical basis ofW'^, the fundamental horizontal vector fields 
are then defined by 

Hi = . 

Now, for every M € Od(K) (space of d x d skew symmetric matrices), we can define a vertical 
vector field Vm on O (M) by 

(yMF)(n) = lim- ^ > ^ ' 



t^o t 

where u e O (M) and F : O (M) -^R. If Eij, I < i < j < d denote the canonical basis of Od(R) 
(Eij is the matrix whose (i,j)-th entry is 1/2, {j,i)-th entry is —1/2 and all other entries are 
zero), then the fundamental vertical vector fields are given by 

It can be shown that we have the following Lie bracket relations: 

k<l 

[Hi, Vjk] = —6ij-Hk + 6ik-Hj, 

where 6ij = 1 if i = j and otherwise, and where is the Riemannian curvature form: 

Q{X, Y){u) = u-^R{tt^X, 7r^Y)u, X,Y € T^O (M) , 

R denoting the Riemannian curvature tensor on M and vr the canonical projection O (M) M. 
In this setting, the Bochner 's horizontal Laplace operator is by definition the operator on O (M) 
given by 

d 

^o{M) = y^.Hf. 

i=l 

Its fundamental property is that it is the lift of the Laplace- Beltrami operator Am of M. That 
is, for every smooth / : M ^ M, 

^0(M)(/ o vr) = (Am/) o vr. 

For the reader unfamiliar with the above construction, we refer for instance to Chapter 3 in 
for further details. It is also proved in the last reference that if the curvature form J7 is 
everywhere non degenerate, then A^^m) 'is subelliptic. Under this last assumption, it is then 
readily checked that the study of A^^m) /fl^^s into our framework. 
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Example 2.4 [The subelliptic Laplace operator on CR manifolds] Let Wi he a non degenerate 
CR manifold of real hypersurface type and dimension d + 1, where d = 2n. Let 9 he a contact 
form on M with respect to which the Levi form Lg is positive definite. Let us assume that the 
pseudo-Hermitian torsion of the Tanaka- Webster connection of (M, 9) is zero. We denote by 
T the characteristic direction of 9 and consider a local orthonormal frame Ti,...,Tn, that is 
Lg{Ti,Tj) = eij where eij is the Kronecker symbol. The following commutations properties hold: 

n 

m,T,] = Y,{r':,-r%)n, 

k=l 

n 

[T,,Tj] = -2V^e,,T + ^ r|i^ - rf^Ifc, 

k=l 

n 

fc=i 

where the T^^ are the Christoffel symbols of the Tanaka- Webster connection (see pp. 32). 
The sub-Laplacian (locally) reads: 

n 
k=l 

where T* is the adjoint ofTi with respect to the volume form 9 A {d9)^ . If we denote 

Xi = -^{Ti + fi),Xi+n = ^^(Ti -f^i). l<i< n, 
then we fall into the previous framework. 
2.2 Canonical connection 

On M there is a canonical connection associated with the geometry of the differential system 
generated by the vector fields Xi, X^. We consider the degenerate metric tensor g on M, such 
that {Xi{x), ...,Xfi{x)} is orthonormal at each point x G M, and for which the spaces 7i{x) and 
V(x) are orthogonal, and g/v{x) = 0- 

Proposition 2.5 On M, there is a unique affine connection V that satisfies the following prop- 
erties: 

• V(7 = 0; 

• V XiXj is horizontal for 1 < i, j < d; 

• \1 Zmn = 0, 1 < m,n < d; 

• If X,Y are horizontal vector fields, the torsion field T{X, Y) = — VyX — [X, Y] is 
vertical and T{Xi, Zmn) = 0, l<i<(i, l<m, n<f). 
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This connection is characterized by the formulas: 

d 



Vx.X, = ^ r^_,.X, = ^ i (4 + u;l - u;},) X,, (2.12) 



k=l k=l 

d 



Vz„„Xi = - ^ ^imn^l^ (2-13) 



£=1 

VZ„„ = 0, (2.14) 
where we have denoted by T^j the Christoffel symbols of the connection. 

Proof. First, it is easy to check that the connection given by the formulas (I2.12p . ()2.13p . ()2.14p 
satisfies the above properties. 

Let now V be an affine connection that satisfies these properties. One first has 



k=l 



Next, by using the fact that = and that the torsion tensor has to be vertical, we easily 
obtain that the Koszul identity holds for V, that is 

g{VxY,Z) = ]^{Xg{Y,Z) + Yg{Z,X) - Zg{X,Y) (2.15) 
+ gi[X,Y],Z) - gi[Y,Z],X) + gi[Z,X],Y)}. 

If we use the orthonormahty assumptions on the X'-s and the Z^nS, and (j2.3p . we obtain 

r'^^=g{Vx,Xj,Xk) = ^{g{[X„Xj],Xk) - g{[Xj,Xk],X,) + g{[Xk,Xi],X,)} 
1 



+ -^]k\ ■ 



Similarly, since T{X, Zmn) = and V XiZmn = one has 

d 

Vz_X, = ^5(Vz_X„X,)X,. 

Using (j2.15p again, along with (j2.4p and (|2.5p . we find 

5(Vz_X„ X,) = ^ {g{[Z„,n, X,IX,) - g{[X„Xe],Z„,n) + g{[Xe, Zmn],Xi)} 



- I —^^ -I- A* \ = —^^ 



zmn • 



□ 

In the sequel it will be useful to have the expression of the torsion tensor on the vector fields 
Xi, Xd. 
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Proposition 2.6 For every k,l = 1, ...,d one has 

m,n=l 

Proof. One has 

d t) 

s=l m,n=l 

Using (j2.12p and the skew-symmetry of the matrix [<-i-'|j]jj=i,...,d we now easily conclude that 

^Ik ~ ^k£ ~ ^£k = 0- 

n 

We also record the following consequence of (|2.12|) 

d 

Vx,X, = -Y,^hXj. (2.16) 

Remark 2.7 We can observe that 



L = Y,Xf-Vx,X, 



i=l 



Example 2.8 For the Example \ 2. 11 it is readily checked that V is the Levi-Civita connection. 

Example 2.9 Let us now identify V for the Example \2.3[ whose notations are in force in what 
follows. Let us recall (see for instance Chapter 3 in f^) that the Ehresmann connection form a 
on O (M) is the unique skew- symmetric matrix a of one forms on O (M) such that: 

1. a{X) = if and only if X e nO{m); 

2. V^^x) =X if and only if X £ VO{M), 

where 7iO{M.) denotes the horizontal bundle and VO(M.) the vertical bundle. It is then easily 



checked that for a vector field Y on 0{ 

d 



k=l 



Let us observe for later use that if X, Y are smooth horizontal vector fields then we have for the 
torsion: 

T{X,Y) = -Vnix,Yy 



Example 2.10 For the Example \ 2.4\ it is an immediate consequence of Theorem 1.3 in JE 
that V is the Webster- Tanaka connection. 
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2.3 Sub-Riemannian distance and local volume growth 

In sub-Riemannian geometry the Riemannian distance dn of M is most of the times confined to 
the background (see in this regard the discussion in section 0.1 of Gromov's Carnot-Caratheodory 
spaces seen from within in [IS]). There is another distance on M, that was introduced by 
Caratheodory in his seminal paper [16j, which plays a central role. A piecewise curve 
7 : [0, T] ^ M is called subunitary at x if for every ^ E T|i,.M one has 

d 
i=l 

We define the subunit length of 7 as Isij) = T. If we indicate with S{x, y) the family of subunit 
curves such that 7(0) = x and 7(T) = y, then thanks to the fundamental accessibility theorem 
of Chow-Rashevsky the connectedness of M implies that S{x, y) ^ for every x,y £ M, see 
|20j . |54j . This allows to define the sub-Riemannian distance on M as follows 

d{x,y) = inf{4(7) | 7 G S{x,y)}. 

We refer the reader to the cited contribution of Gromov to [lOJ, and to the opening article by 
Bellaiche in the same volume. 

We next recall that a metric space {S, p) is called a length-space, or intrinsic, if for any x,y £ S 

p{x,y) = inf£(7a:y), 

where the infimum is taken over all continuous, rectifiable curves jxy, joining x to y. For a 
continuous curve 7 : [a,b] — > 5, one defines £(7) = supYl^=i p{7{'ti),'y{'ti+i)), the supremum 
being taken on all finite partitions a = ti < t2 < ■ ■ ■ < tp < tp^i = 6 of the interval [a, b]. Since 
by the triangle inequality we trivially have p{x,y) < ii'jxy) for any continuous curve joining x 
to y, it follows that in any metric space {S,p), 

p{x, y) < inf l{yxy), x,y £ S. 

In particular, such inequality is therefore valid in the space (M, d). Now it is proved in Proposi- 
tion 2.2 in [22] that in any sub-Riemannian space the opposite inequality is also valid. Therefore, 
every sub-Riemannian space (M, d) is a length-space. 

Another elementary consequence of the Chow-Rashevsky theorem is that i : (M, d) ^ (M, dji) 
is continuous. On the other hand, it was proved in [52j that for any connected set C M which 
is bounded in the distance d^ there exist C = C{i}) > 0, and e = e(r2) > 0, such that 

d{x, y) < CdR{x, yf, x,y eQ. 

This implies that also the inclusion i : (M, d^) ^ (M, d) is continuous, and thus, the topologies 
of dfi and d coincide. In particular, compact sets are the same in either topology. The metric 
space (M, d) is locally compact and, furthermore, for every compact set K C Wl there exists 
ro{K) > such that for x £ K and < r < tq, the closed balls B{x,r) in the metric d are 
compact, see Proposition 1.1 in [31j. If M is unbounded and the vector fields grow too fast at 
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infinity, then balls of large radii may not be bounded in da in general, and the space (M, d) may 
fail to be complete. The lack of completeness and the fact that not all metric balls are bounded 
in the metric dpi are two equivalent properties in view of the following powerful generalization 
of the classical Theorem of Hopf-Rinow due to Cohn-Vossen, see |12j. 

Theorem 2.11 (of Hopf-Rinow type) In any locally compact length-space {S,p), the com- 
pleteness of the metric space is equivalent to the compactness of the closed balls. 

Since as we have mentioned (M, d) is a locally compact length-space, if we want to guarantee the 
compactness of closed balls (in the d metric) of arbitrary radii, we need to assume that (M, d) 
is a complete metric space. By what has been said, this is equivalent to requiring that the space 
(M, dji) is complete. 

In section [12] it will be expedient to work with yet another distance on M. For every x,y € M, 
we set 

p{x,y)= sup |/(x)-/(y)|, (2.17) 

/ec°°(M),r<x,(/)<i 

where roo(/) = sup r(/, /), and the quantity T{f,f) is defined in (|4.26p below. Thanks to 

M 

Lemma 5.43 in [ITj we know that 

d{x,y) = p{x,y), x,yeM, 
hence we can work indifferently with either one of the distances d or p. 

The following fundamental result proved by Nagel, Stein and Wainger in [52] provides a uniform 
local control of the growth of the metric balls in (M, d). 

Theorem 2.12 For any x G M there exist constants C{x),R{x) > such that with Q{x) = 
log2 C{x) one has 

li{B{x,tr)) > C{x)-^t^^''^p{B{x,r)), 0<t<l, 0<r< R{x). 

Given any compact set K dWi one has 

inf Cix) > 0, inf R(x) > 0. 

2.4 Geodesies 

In this section we prove a first variation formula for the geodesies of the metric space (M, d). 
More precisely, we will describe with our connection V the local minima of the energy functional. 
To this end, let us first introduce a family of natural skew-symmetric linear operators of the 
horizontal bundle. For V £ 7i, we define: 

d 

JmniV) = g{X,, V)i^^Xi, 1 < m, n < f). 

With this in hands, we can now extend the method of Rumin [59] to provide the equation of 
the geodesies. We first have the following result. 
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Lemma 2.13 For 1 < m,n < I), let us denote by 6mn the one-form on M such that 0ran{Zmn) = 
1 ; Qmn{Zpq) = if {m,n} 7^ {p,q} and 0{Xi) = 0. IfVi is a smooth vector field on M, and V2 
is a smooth horizontal vector field, we have 

demn{Vl,V2) = -g{Vi,J^nV2). 

Proof. Let us first assume that Vi and V2 are both horizontal. From Cartan's formula 
Therefore 

de^n{Vl,V2) = -9„,n{[Vl,V2]). 

^From the definition of we easily check that for 1 <i,j < d, 

9mn{[Xi,Xj]) = J^"-, 

SO that 

d9mn{Vl,V2) = -g{Vi,JmnV2)- 

Now, if Vi = Zpq, since from our assumption (j2.4p the bracket [Zpq,Xi] is always a horizontal 
vector field, again from Cartan's formula we obtain 

d9^niVuV2) = 0. 

n 

Proposition 2.14 Let 2;,y G M and denote by T-L{x,y) the set of horizontal curves j{t), < 
t < 1, going from x to y. Then 'j & Ti.{x,y) is a critical point of the energy functional 

E{l)= (\{^'{t),i{t))dt 
Jo 

if and only if there exist constants amn, 1 < m, n < [) such that 

m,n=l 

Proof. Let us consider a family of curve 7u(t), —€< u < e, < t < 1, such that 7o(t) S Ti.{x,y) 
and V = is a Legendre vector field of the constrained problem corresponding to the 

variational problem we are looking at. We denote X = j'(t). We have 

V{x) = V{y) = 0, 

and 

[V, X] G n. 
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Since the torsion of V is vertical, we now compute 



= Cvg{X,X)dt 
du |«=o Jo 

= 2 / g{VvX,X)dt 
Jo 

= 2 / giVxV,X)dt 
Jo 

= 2 / Xg{V,X)-g{V,VxX)dt 
Jo 

= -2 [ g{V,VxX)dt 
Jo 



Therefore, if 7 is a critical point, we must have g{X, VxX) = 0. We now write 

t) 

^xX = ^ amnit)JmnX + Y, 
m,n=l 

where Y is orthogonal to X and to Ylmn=i '^mn{t)Jmn^ ■ We have then 

{t)g{V,JmnX). 

m,n=l m,n=l 

But, from the previous lemma and Cartan's formula 

g{V, JmnX) = -dOmniV, X) = Xe{V). 

Therefore, 

^^P^ =-2 /' E amn{t)Xe{V)+g{V,Y)dt 

(^U \u=0 Jo 

Integrating by parts, we obtain 



-2 / X -a'mn{mV)+9{V,Y)dt, 
•^0 m,n=l 



du \u=0 

so that 7o is a critical point if and only if is constant and Y = Q. □ 
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3 The curvature tensor 



In this section we introduce a first-order differential quadratic form wliich plays a pervasive role 
in tlie results of tliis paper. If / G C°°(M.) we define 

d f / d i) \ d 

k,e=l ^ ^ j=l m,n.=l ^ j=l 

jj=l i=l l<i<j<d ^ ^ 

+ E E ( E + E - E ^.-^r) 

fc=lm,?i=l ^£,j=l l<i<j<d j=l ^ 

\ 2 



+ ^ E ( E 7r^--/) • 



\<l<j 

The geometric meaning of the differential quadratic form /) will become fully clear in 
section HI The main purpose of the present section is to prove its tensorial nature. In fact, in 
Proposition 13.11 below we show that '7^(/, /) can be expressed solely in terms of the curvature 
and torsion tensors with respect to the canonical connection V introduced in section [2^21 At the 
end of this section we work out the formulas for /) in some special examples. In example 
13.21 we show that in the Riemannian case this form coincides with the Riemannian Ricci tensor. 
In what follows V indicates the canonical connection on M introduced in section 12.21 Let us 
recall that the torsion tensor is given by 

T(X, Y) = VxY - VyX - [X, y], 

the curvature tensor is given by 

R{X, Y)Z = Vx^yZ - Vy^xZ - Vyx,Y\Z, 

and the Ricci tensor by 

d 

Ric(X,y) = Yg{R{Xi,X)Y,Xi), 

i=l 

where X, Y, Z are smooth vector fields. 

Proposition 3.1 If f : M. ^ M. is a smooth function, then 
d 



T^UJ) = Yl ^HXi,Xk)XefXkf - {{Vx,T)iXe,Xk)f)iXkf) + - {J{X,,X^)fY 



i,k=l 

As a consequence, IZ is a (0, 2) tensor on M and therefore it is coordinate free. 
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Proof. We begin by writing the quantity TZ{f, f) in (|3.18p as follows 

nf, f) = T^iif, f) + T^iiif, f) + T^iiiif, f), 

where 

k,£=l ^^j=lm,n=l ^ j=l 

i,j=l i=l l<i<j<d ^ ^ 

T^Ilif, /) = E E ( E + E 47.7" - E ^^-7.7) ZmnfXkf, 

k=lm,n=l ^e,j=l l<i<j<d i=l ^ 



T^iii{fJ) = \ E ( E 7^^)' 

l<l<j<d ^m,n=l ^ 



l<^<i<d 

The proof will be completed if we show that 

d 

T^iUJ) = E Ric(Xfc,X,)Xfc/X,/, (3.19) 

d 

nnifj) = -Y {{^x,T){X,,Xk)f){Xkf), (3.20) 
e,k=i 

d 

/) = 7 E (^(^^' ^'^)/)' • (3.21) 
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To establish (I3.19P we observe that 



d d t) d 

mc{Xk,Xe) = E 9iRiX^, Xk)Xe, ^.) = E E ^h^^Un + E {^:>^ie " ^'^r^'.) (^.22) 

j=l j=l m,n=l j=l 

d 

The identity (I3.22p can be verified in a standard fashion. Our goal is thus proving that 

E ( E (^ri, - XkT^ + Y (ri,n, - n^l, - u^lk^j?! ] {Xkf){x,f) 

k,£=l \j=l i,j=l J 

d / d d d 



a / a a a 

E ( E(^^<- - ^^-4) + E ^ii^fci - Y'^ki^ 

k,i=l ^ j=l i,j=l i=l 

+ 1 Y - K- + + O) ) {Xkf){Xef). 



\<i<i<d 
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To show this we use the formula (|2.12p 

= - ( ^^ij + - J ) 

which allows to express the Christoffel symbols F^^- in terms of the structural constants cufj. 
Keeping in mind that cofj = —^^ji, one easily recognizes that 

d / d 

k,i=i \j=i ) k,i=i \j=i 

To complete the proof of (j3.19p we are thus left with proving that 

d d 

E E i^ii^h - r^r^, - u;^rj,) Xkfxj (3.23) 

k,e=i i,j=i 

= E { E^H.-E'^fc«^^^ + ^ E 

k,e=l ^ i,j=l 1=1 l<i<j<d ^ ^ 

We now notice that again from (j2.12p and the skew-symmetry of wj^ in the lower indices we 
obtain 

E E n^h^^m = E E ^^^n,Xkfx,f 

k,e=l i,j=l k,e=l i,j=l 

c( d ^A: _|_ d d 

= E E '\ '' ^h^kfxd = E E ^%^\,Xkfxd 

k,£=li,j=l k/=li,j=l 
d d 

= E E ^ir^jXkfxd. 

k,i=l i,j=l 

Next, we have 

d d d d 

-EE (r^r^, + ^^q,) Xkfxd = - ^ ^ (r^^ + u;j,) Xkfxj. 

k,i=l i,j = l k,i=l i,j=l 



By ([21^ one has 



rij + ^ik = 2 y^lk + ^ij + ^jk 
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and thus 

d d 



:/=l jj=l 
d d 

7 E E (^^ - - ^h) (4 +^lk+ xufxd 



4 

k,f.=l i,j=l 



d d 1 
k,e=l 1=1 k,e=l l<i<j<d 

IE E {^h+<)H+^L)XkfXef, 



k,i=l l<i<j<d 

where the last equahty is obtained by expanding the product and cancehng equal terms. This 
proves ()3.23p . thus completing the proof of ()3.19p . 
Next, we turn to the proof of ()3.20p . One has 

d d 

{iVx,T)iXe,Xk)f)Xkf = Vx, {T{X,,Xk)) fXkf (3.24) 

i,k=l l,k=l 
d d 

- E T{Vx,Xe,Xk)fXkf- J2 T{Xe,Vx,Xk)fXkf. 

i,k=l e,k=i 

Proposition 12.61 and the fact that VZ = now give 

d d i) d d d 

Y VxAnXi,Xk))fXkf = -Y E {XaT)ZmnXkf = Y E E(^^-^D^-"^fc/- 

£,k=l £,k=l m,n=l k=l m,n=l j=l 

Next, using (|2.16p and Proposition 12.61 again one obtains 

d d d 

- E T{Vx,X,,Xk)fXkf= Y Y.'^ljT{Xj,Xk)fXkf 

e,k=i e,k=ij=i 

= -E E E47r^-/^^/- 

k=l m,n=l j,i=l 
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Finally, 

d 



T (X,, Vx.Xk) fXkf = - Yl ^ikT (Xe, Xj) fXkf 

t,k=l j,t,k=l 

E E t.^T^lkZn.nfX,f = lY. E Y.^T(^lk + ^'i + ^%)ZmnfX,f 



2 

k=l m,n=l j ,£=1 k=l m,,n=l j/=l 

d t) 

2 



^ E E E ^T(^k + 4 + - - 4 - OZmnfXJ 
k=l m,n=l l<l<j<d 
d fl 

E E E 7r4^-"/^^/- 

A;=l m,n=l l<l<j<d 



Substitution in (j3.24p gives 

d 



{{Vx,T){Xe,Xk)f){X,f) = -TZnifJ), 



£,k=l 



which proves (13.201) . In order to complete the proof we are left with establishing (I3.2ip . From 
Proposition 12.61 and the skew-symmetry of 7^" in the lower indices, we easily find 

d d / t, Y 

-Yinx,,x,)ff = -Y E =njn{f,f). 

e,k=l e,k=l \m,n=l J 



This finishes the proof. 

We can now compute TZ in the geometric examples introduced in section \27L[ 
Example 3.2 For a Riemannian manifold M, see Example \2.1[ we have 

7^(/,/)=Ric(V/,V/). 
This follows immediately from Proposition \3.1[ 



□ 



Example 3.3 When M is a graded nilpotent Lie group of step two, see Example \2.^ we imme- 
diately obtain from (j3.18p 



W,/) = \iZ^z,,ff = - Ym.x,]f)\ 



ij=i «j=i 

We note here that, in terms of the quadratic form introduced in (j4.29p below, one has 

7^(/,/) = V(/,/). (3.25) 
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Example 3.4 In the Example \2.3\ we have 

12 

7^(/,/) = ^ mc*{Hj,H,){Hjf){Hkf) + Vy„^n(H„H,)fHkf + -{Vn^H„H,)f) , 

j,k=i 

where for horizontal vector fields X and Y , 

Ric*(X,y) = Ric(^,X,7r,y), 
with Ric denoting the Ricci tensor o/M. 

4 Bochner type formulas 

Our goal in this section is to prove sub-Riemannian Bochner type formulas which play a central 
role in this paper. Since the operator L is two-step generating, we prove two types of such 
formulas, one involving the commutation between the horizontal gradient and L, the other one 
involving the commutation between a vertical gradient and L. As it is to be surmised, the 
horizontal one will prove quite involved technically. 

In the sequel, we will use the following differential bilinear forms: 

1 ^ 
T{f,g) = -{L{fg) - fLg - gLf) =Y,XifXig, (4.26) 

i=l 

r2(/, g) = \ [LT{f, g) - T{f, Lg) - T{g, Lf)] , (4.27) 
T^{f,g)=Y,iZ,,f)iZijg), (4.28) 

rf (/, g) = ^ [LT^U, 9) - r^(/, Lg) - V\g, Lf)] . (4.29) 

4.1 The first Bochner formula 

Henceforth, we adopt the notation 

f _ XjXjf + XjXjf 

J,ij ~ 2 

for the entries of the symmetrized Hessian of / with respect to the vector fields Xi, X^- 
Noting that 

XiXjf = f^ij + - [Xi,Xj\f, 

using (j2.3p we obtain the useful formula 

. d t) 

XiXjf = + + - ^ iT^rnnf ■ (4.30) 

£=1 m,n=l 

Our principal result of this section is the following: 
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Theorem 4.1 (Horizontal Bochner formula) For every smooth function / : M — > M, 

T2{fJ)=j:ifM-Y.-'^X,f] +2 E f - E ^^^^^^«/) (4-31) 

1=1 \ i=l / l<«j<d \ i=l / 

d f) 

-^E E lTi^3Zmnf){Xd)+n{fJ), 
i,j=l m,n=l 

where TZ{f,f) is the quadratic form defined in (13.18^ . 
Proof. We begin by observing that for any smooth function F on M 

= 2FLF + 2T{F,F), 

This and (OGI) give 

d d d 

LT{f, /) = E L{{Xiff) = 2 X.JLiXif) + 2 nx^f, X,f). 

i=l 1=1 i=l 

We now have 

d d 

L{X,f) = XoXif + Y X]XJ = XiXof + [Xo,Xi]f + ^ Xj{XiX,f) + Xj[X,,X,]f 
i=i i=i 

d 

= Xaof + [Xo,X,]f + Y {X^{X,X,f) + [X,,X,]Xjf + X,[X,,X,]f} 

= X,{Lf) + [Xo, X,]f + Y { [Xj,Xi]Xjf + Xj[X,,X,]f} 

i=i 

d d 

= X,{Lf) + [Xo, X,]f + 2Y,[Xj,Xi]X,f + [Xj,Xi]]f. 

Using this identity we find 

d ( d d 

LT{f, f) = 2Y,Xifl X,{Lf) + [Xo, X,]f + 2Y,[Xj,Xi\Xjf + ^^[X,-, [X,-, X,]]/ 
i=i [ j=i j=i 

d 

+ 2 (X,XJ)2 

d d d 

= 2r(/, L/) + 2 XJ[Xo, X,]/ + 4 X,/[X„ X,]X,/ + 2 ^ X,/[Xj, [X^, X,]]/ 

i=l i,j=l i,j=l 

d 

+ 2 Y (XjXiff. 
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Since, thanks to the skew-symmetry of the matrix {[^j, Xj]/}ij=i,...,d, we have 

d 

Y^f,,[X„X,]f = 0, 

we find 

d 



4 



We thus obtain 

- [Lr(/, /) - 2r(/, L/)] = x^f[Xo, x,]f 

1=1 

d d 

+ 2 Xif[X„Xi\X,f + J] 

i,j=l i,j=l 
ij=l jj=l 

Since 623) gives Tal/, /) = i [Lr(/, /) - 2r(/, L/)] , we conclude 

d d 

r2(/, /) = E - 2 E XJ[X„X,]Xjf (4.32) 

ij=l jj=l 



^ a a a 

To complete the proof we need to recognize that the right-hand side in (j4.32p coincides with 
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that in (j4.3ip . With this objective in mind, using (|2.3p and ()4.30p we obtain 
d d 

i,j=l i,j=l 

d d f d i) 

1=1 l<i<j<d i,j=l \l=l m,n=l 

= E/,^« + 2 4. 

£=l l<i<j<d 

d d d i) 

- 2 E E - 2 E E ^T^rnnXjf XJ 

i,j=l i=l i,j=i m,n=l 

d d d 

= E/,« + 2 E /'^-2EE4/a-^^/ 

e=i i<i<j<d i,j=i e=i 

-EE ^f.^-^'^/^^/ - E E E -^f^^T^rnnf XJ 
i,j=l £,k=l * ,i=l ^=1 m,n=l 

d i) 

- 2 E E ^ij'^^mnXjf Xif 

i,j=l m,n=l 

= E/,l. + 2 E /,?.-2E fE4-x./)//. 

£=1 i<t<i<d e,j=i \i=i / 

-EE - E E E ^f.7r^-/ xj 

i,j=l t,k=l i,j=^ ^=1 m,n=l 
d i) 

- 2 E E ^ij^^rnnXjf Xif 

i,j=l m,n=l 



^<^<j<d \ l<^<j<d \i=l / 

-EE -^fA^>^f^^f - E E E ^f.Tr^-'^/ 

i,j=l£,A:=l «,J=1 ^=1 m,n=l 

- 2 E E ^ij'^^rnnXjf Xif. 
i,j=l m,n=l 
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If we now complete the squares we obtain 

d d 

f% - 2 XJ[X„X,]X,f (4.33) 



i,j=l i,j=l 

1=1 \ 1=1 / l<f.<j<d \ i=l 



2 



£=1 \i=l / l<l<j<d \i=l , 

d d d \) 

- E EE E ^pi^'^z^nf 

i,j,k,£=l « J=l ^=1 m,n=l 

d h d h 

-2E E iT^jZrnnf Xif-2Y, E 7."" [^-- ^^■]/ ^^Z" 
j j'=l m,n=l *ij=l m,n=l 

Next, we have from (|2.7p 

E^./[^o,^.]/= E (4.34) 

j=l i,j,k,£=l 

d d f) d d 

+ E E E ^^fc7r^--/^^/ + E E (^.^J^)^./^./, 

j=l j,k=l m,n=l i=l j,k=l 

and also 

d d d d i) 

E = E EK-^^'^^-]/^^/ + E E [7r^™«'^^-]/^^/ 

i j=l i J=l i=l i,j=i m,n=l 

d d d d 

= E E4-^^/[^^'^^-]/- E E(^^-^^.?-)^*/^^/ 

ij=i e=i i,j=i 1=1 

d f) d f) 

+ E E ^rt^-"'^^-]/^^/- E E (^^-^D w^^/. 

j j=l m,n=l * >i=l m.,»i=l 

Using (j2.3p we find 



d d d d f) 

E^j[[^.,^,],^,]/= E E^'.4^*/^'^/+ E E 4-7^^--/^^/ (4-35) 

i,j=l «,i=l ^,fc=l jj',£=l m,n=l 

d f) d f) 

+ E E ^rt^--^^-]/^^/- E E (^^-^D w^./ 

i,j=l m,n=l i,j=l m,n=l 

d d 

- E E(^j-^^j-)^^/^^/- 

i,j=l £=1 
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Again by (|2.3p we have 



^i:([^.x,]/)^ = ^ E (E^^^f) +\ ^ fi^rw) (4.36) 

i,j=l l<i<j<d \l=l J l<i<j<d \Tn,n=l J 



l<i<j<d e=l m,n=l 



Substituting (|i33|) - (06]l in (021) we obtain 

r.(/,/) = Ef/,«-E-f.^^/) +2 E f/.^-E^^^^^/V 

£=1 \ 1=1 / l<t<j<d \ i=l / 

d f) 

-2 E E iT ^ J f Xa + Monster 

i,j=l m,n=l 



where we have let 



d / d \2 / d e , j \2 

Monster = -Y^ I Y^u^'^Xj] -2 E E (4-37) 

£=1 \j=l / l<i<j<d \i=l / 

d d d d t) 

+ E ^^Ujf^XefXd- E u;^A,^efXd-Y.Y.J2^'j^T^rnnfXif 

i,j,k,£=l i,j,k,£=l i,j=i (-=1 m,n=l 

d tf d d t} 

-EE 7r[^--^^]/^^/+E E E ^jkir^-^nfxj 

i,j=l m,n=l i=l j,k=l m,n=l 

d d d d d i) 

+E E(^*^''fc)^^/^j-/+ E E4-4^^/^^-/+ E E 

4=1 j,fc=l i,j=ll,k=l i,j,i=lm,n=l 

d t) d d 

-EE (^,7r)^"-/^^/- E E(^^-4)^«/^^/ 

i,j=l m,n=l 'i'd=^ £=1 

2 



+ \ E ('e4-w)%^ E (e 

l<i<j<d \£=1 / l<i<j<d \m,n=l 

+ E E E ^hTZmnfXJ. 

i<i<j<d £=1 m,n=l 



mn y f 
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Simplifying the expression we obtain 

d d 



Monster = - ^ ^ w^^w^^Xfc/X^/ 

k,l=l i=l 

1 

-oE E (4.38) 



2 

fc,;=l l<i<j<d 

d d d 

k,e=lj=l i,j,k,i=l 
d d i) 

+ 2E E 44^^-w+E E ^rt^^'^^"]/^^/ 

fc/=l l<i<j<d k,j=l m,n=l 

+ E E E ^^;.7.'^"^^n/X./ + E E E -^fjlT^rnnfXd 
2=1 j,k=l m,n=l ^<i<j<d ^=1 ni,n=l 

d f) 1 / " \^ 

- E E (^^■7r)wxj+- . 

i j=l m,?i=l l<i<j<rf \m,n=l j 

To complete the proof we need to recognize that the right-hand side of (j4.38p coincides with 
7^(/, /) defined by (f3J8]) . If we now use ([23]) we obtain 

d fl d ^ d 

E E 7feT[^„^™n]/^fc/= E E^fcT'^l-n^'^/^^/ 

k,j=lm,n=l k,j=l'm,n=l £=1 

d f) 



E (E E T.T'^in.n)^^./^./, 
A:,£=l ^ j=l m,n=l ^ 



Substituting in (j4.38p we find 



Monster = E |( E E 7.7^'-,) + E(^^<- " X,uj%) (4.39) 
+ ^^,4^.-^u;^,4 + i ^ (44-(4,.+a;y(4^.+4)]|x,/XJ 



i,j=l i=l l<'<j<ci 

+ E E E ^)kiT^^nfxd + E EE 4^rzmnfxd 

1=1 j,k=l m,n=l l<i<j<d. 1=1 m,n=l 

d i) 1 / ^ \^ 

-EE iXjlD^mnfXJ + 2 E E ^r^rnnf ■ 

i,j=lm,n=l l<i<j<d \m,n=l j 
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Rearranging the indices we can rewrite (j4.39p as follows 

d f / d i) \ d 

Monster = E E E 7^7^'- + E(^^<- " ^^■4) 

kl=l ^ ^ 1=1 m,n=l ^ j=l 



d 



ij=l i=l l<i<j<rf 

+EE(E47.t+ E 47r-E^.7r)w^./ 

+ ^ E ( E iT^rnnfX =n{fj). 
l<e<j<d ^m,n=l ^ 

This completes the proof of Theorem I4.1i □ 

Remark 4.2 When M. is a graded nilpotent Lie group of step two, then using the structural 
constants in Example \2.S\ we obtain from formula ()4.3ip 

d d f) 

r2(/,/) = E/,« + 2 E f%-^^ ^ iT(^jZ^nf){xd) + nf,f). 

e=l i<^<j<d i,j=lm,n=l 



Keeping ()3.25p in mind, we thus we obtain from the latter equation 

d d 

4 



r2(/,/) = ||v|,/||2 + - E {[x^,x,]ff + 2 E 



i,j=l i,j=l 

which is Proposition 3.3 in JSU] (the latter, however, is relative to groups of arbitrary step). 
Here, V|^/ = [f^ij] represents the symmetrized horizontal Hessian of f . 

4.2 The second Bochner formula 

To prove our next result we will need the following lemma. 

Lemma 4.3 For every m,n = 1, f) one has 

[L, Z„in] = 0. 

Proof. We begin by observing that, as consequence of the assumption (j2.5p . for every m,n = 
1, t), the commutator[L, Zmn] is a vector field. Indeed, given a smooth function /, using the 
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formulas ()2.4p . ()2.7p and ()4.30p . one obtains by elementary computations 

d 

[L, Zmn]f = [Xq, Zmn]f + Xi[Xi, Zmn]f + [Xi, Zmn]Xif 

i=l 

d d d 

i,k,i=l i,k=l i,l=l 

d d d 

i,fc/=l i,fc=l i,£=l 



where we used the crucial fact that 
To complete the proof let 4i,tp £ C^(M), then (p^ gives 



< [L,Zm„]>,'0 > =< 0, [L,Zmn]V' > = < 0, -^^(^mnV') > " < (p, ZmniLtp) > 

= < L(p,Zmn1p > + < Zmn(i),L'4) > = < L{Zmn4>,'4' > " < ZmniL(p),'4> > 
= < [L, Zmn]4>,'^ >, 

and thus [L, Z^n]* = [L, Zmn], i-e. ^mn] is a symmetric vector field. Such vector fields must 
vanish. Indeed, if V is symmetric vector field, one has for any (jj^ip & C°°(M) 

= < >=< V*IA^ >=< ^,V{(l)^p) >=< l,(j)Vij > + < l^^jVcp > 

=<(/), > + < V', V(i)>=2 <V(i),ip > . 

Taking ^ = Vcj) we conclude that it must be Vcj) = 0. By the arbitrariness of (/) G C°°(M) we 
conclude ^ = 0. 

□ 

Our second Bochner type formula is expressed by the following proposition. 
Proposition 4.4 (Vertical Bochner formula) For every smooth function / : M ^ M, 

rf(/, /)= T(Zmnf, Zmnf) 

m,n=l 

Proof. ^From (j4.29p we obtain 

1 



rf(/,/) 2 



LT^{fJ)-2V^{f,Lf) 
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We now have from (j4.28p and from Lemma 14.31 

t) f) I) d 

LT\fJ)= L{{Zmnff) = 2 J2 Zmnf L{Zmnf) + 2 

m,n=l m,n=l m,n=l k=l 

f) f) d 

= 2 Zmnf Zmn{Lf) + 2 (^Xk{Zmnf)) 

m,n=l m,n=l k=l 

= 2r^(/,L/) + 2 ^ r(z„„/,z„„/), 

m,n=l 

where in the last term we have used (|4.26|) . □ 

5 A sub-Riemannian curvature dimension inequality 

As a consequence of our two Bochner's formulas we obtain a generalization of the curvature 
dimension inequality of Bakry. A new feature of such inequality is the presence in the left hand- 
side of the vertical quadratic form Ff (/, /). For a smooth function /, we define the quantity 

d t) / d \ ^ 

j=l m,n=l \«=1 / 

We observe that with the notations of Section [2 



T{f,f)= E 5(^n^n(V/),Jn^„(V/)). 
m,n=l 

We notice explicitly that, since in the Riemannian case 7^" = 0, see Example 12.11 in that case 
we have 

T{fj) = 0. (5.41) 

When instead M is a graded nilpotent Lie group of step two, then from the expressions of the 
7ij " = in example 12.21 we find 

m/) = ^r(/,/). (5.42) 

Proposition 5.1 (Sub-Riemannian curvature-dimension inequality) For every smooth 
function / : M ^ M and every > 0, 

F2(/, /) + i^Tlif, f) > -(Lff + n{f, f) - -Tif, /). 

a V 

Proof. /,From ()2.6p and ()2.7I) and Schwarz inequality we find 

d / d \ ( f X 2\ 1/2 

= E - E^^^^^/ U ^ E -Y.^lXif 

1=1 \ 1=1 J \^=l \ i=l / 
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^From this inequality and from Theorem 14. H using Schwarz inequality again we obtain for every 

d h / d \ 

7=1 m,n=l \i=l / 



2 E 



j=lm,n=l j=lm,n=l \i=l / 

Using Proposition 14.41 and the definition (j5.40p we find 

^{Lff < T^if, f) + i^Tlif, /) + ir(/, /) - 7^(/, /), 
which gives the desired conclusion. □ 

Example 5.2 On a Riemannian manifold M with Laplace- Beltrami operator L we have (/, /, ) 
T{f,f) = (see (fOT]) ). andn{f,f) = Ric (V/, V/) (see ExamplelEE). PropositionlM thus 
gives 

T2if,f)>^{Lff + mc (V/,V/). 

When 

Ric{V,V) > pi e M, (5.43) 

i(;e thus recover the Riemannian curvature- dimension inequality 

r2ifJ)>^{Lff + p^T{fJ). (5.44) 



Remark 5.3 It is important to keep in mind that, in view of Theorem 1.3 in |57| / and Proposi- 
tion 3.3 in f4^, the Riemannian curvature dimension inequality (j5.44p is, in fact, equivalent to 

Example 5.4 Let us assume that L is the sub-Laplacian on a two-step graded nilpotent Lie 
group, then for every smooth function f and every v > 0, 

r^if, /) + ^rf (/, /) > i(L/)2 - ^r(/, /) + ir^(/, /). 
6 Gradient estimates for the heat semigroup 

In this whole section we assume that the metric space (M, d) is complete, where d denotes 
the Caratheodory distance associated to L (see section r2.3p . Recall that this is equivalent to 



33 



assuming that M is complete with respect to its Riemannian distance dn. We also suppose that 
there exist constants pi G M, p2 > and k > such that for every smooth function / : M — > M: 

n{f,f)>piT{fj) + p2T^{fj), 

T{f,f)<KT{f,f). 

We emphasize that, as we have seen in example 15.21 and remark 15.31 the assumptions (j6.45p 
should be considered as the sub-Riemannian analogue of (j5.43p . In particular, the requirement 
Pi > in (j6.45p corresponds to the Riemannian Ric> 0. 

According to Proposition 15.11 the assumptions ()6.45p imply, for every smooth function / : M ^ M 
and every > 0, 

r2(/, /) + z.rf (/, /) > i(L/)2 + (pi - ^) r(/, /) + p2r^(/, /). (6.46) 

This is a sub-Riemannian version of the curvature-dimension inequality (jl.ip . 



6.1 The heat semigroup 

As an operator defined on C^(M) the operator L is symmetric with respect to the measure p 
and non positive: For / G Co"(M), < Lf, / >< 0. Therefore, it admits a self-adjoint extension 
to L^(M, p), the Friedrichs extension. Following an argument of Strichartz [M], Theorem 7.3 p. 
246 and p. 261, we now prove that L is essentially self-adjoint on C^(M) (the completeness of 
the metric space (M, d) is crucial here). 

In what follows, to distinguish it from the canonical connection V introduced in section 12.21 
we use the notation V/j for the Riemannian connection on M. We recall that, thanks to the 
assumption (j2.1ip . 

IIVk/IP = r(/) + V(/), /GC°°(M). 

Lemma 6.1 There exists an increasing sequence hn € C^(M) such that hn 1 on M, and 
||V/j/i„||oo ^0, as n ^ oo. 

Proof. As in [32j, if we fix a base point xq £ M, we can find an exhaustion function p G C°°(M) 
such that 

\p-dR{xo,-)\<L, \Vrp\<L on M. 

By the completeness of (M, dji) and the Hopf-Rinow theorem, the level sets Jig = {x G M | 
p(x) < s} are relatively compact and, furthermore, $7^ M as s — > oo. We now pick an 
increasing sequence of functions cpn G C°^{[0,oo)) such that = 1 on [0,n], = outside 
[0,2n], and < f . If we set /i„(x) = 0„(p(x)), then we have /i„ G C^(M), /i„ /" 1 on M as 
n — > oo, and 

2L 

W^RhnWoo < • 

n 

□ 

In what follows we define the action of the operator —L on C^(M) by the equation 

<-L(t>,ij>= [ T{(f>,'^)d^l, -^eC^iM). 

Jm 

We recall our assumption (|2.9p . As a corollary of the previous lemma, we obtain: 
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Proposition 6.2 The operator L is essentially self-adjoint on C^(M). 

Proof. Denote by L the Friedrichs extension of L initially defined on C^(M). If L* is the adjoint 
of L, then according to Reed-Simon [SB], p. 137, it is enough to prove that the eigenvalues of 
the adjoint L* are negative. We thus have to show that \i L* f = Xf with A > 0, then / = 0. 
From the hypoellipticity of L, we first deduce that / has to be a smooth function. Now, for 
h e C^{M), 

< Tif, h^f) >=-</, >= - < L*f, h^f >= -A < /, h'^f >= -X < f\ >< 0. 

Since 

T{f,h^f) = h'T{fJ) + 2fhT{f,h), 

we deduce that 

</i2,r(/,/)>+2<//i,r(/,/i)><o. 

Therefore, by Schwarz inequality 

<h\T{fj) ><ml\\nKh)\w 

If we now use the sequence /i„ constructed in Lemma l6. II and let n — > oo, we obtain r(/, /) = 
and therefore / = 0, as desired. □ 

If L = — /q^°° XdE\ denotes the spectral decomposition of L in L^(M, /u), then by definition, the 
heat semigroup {Pt)t>o is given by Pt = /q^°° e^'^^dEx- It is a family of bounded operators on 
L^(M, //). Since the quadratic form — < f,Lf > is a Dirichlet form in the sense of Fukushima 
|29j . we deduce from the above result that C^(M) is dense in the domain of Pt and that (Pt)t>o 
is a sub-Markov semigroup: it transforms positive functions into positive functions and satisfies 

Pjl < 1. (6.47) 

This property implies in particular 

\\Ptf\\L^(M) < II/IIli(M): Il-Pi/I|L°°(M) < II/IIl°°(M)> (6.48) 

and therefore by the Theorem of Riesz-Thorin 

II^'JIIlp(m) < II/IIlp(m), 1<P<oo. (6.49) 
Moreover, it can be shown as in [63j, Theorem 3.9: 
Proposition 6.3 The unique solution of the Cauchy problem 

ff -Ln = 0, 

\n(x,0) = /(x), / G LP(M),l<p < +CX), 
that satisfies ||ii(-,t)||p < Ce*^*, for some constants C and M, is given by u{x,t) = Ptf{x). 
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Due to the hypoellipticity of L, (t, x) Ptf{x) is smooth on M x (0, oo) and 

Ptf{x)= [ p{x,y,t)f{y)dfi{y), feC^{M), 



where p{x, y,t) > is the so-cahed heat kernel associated to Pf. Such function is smooth outside 
the diagonal of M x M, and it is symmetric, i.e., 

P{x,y,t) =p{y,x,t). 

By the semi-group property for every x, ?/ G M and < s,t we have 

Pix,y,t+s)= / p{x,z,t)p{z,y,s)dfi{z) = / p{x, z,t)p{y, z, s)dniz) = Ps{p{x, ■,t)){y). (6.50) 
Jm Jm 

6.2 A variational inequality 

The goal of this section is to establish a variational inequality which is the cornerstone of all the 
gradient estimates we shall obtain in the sequel. Henceforth, to simplify the notation we let 

r(/) = r(/,/), fGC\M). 

For the sequel it will be convenient to observe that if G C^(M), / G C^(M), then 

L(c/> o /) = 0"(/)r(/) + cP'if)Lf. (6.51) 
In particular, if / G C^(M) and / > 0, then for any fixed T > and t < T, we obtain from 

LPr tf 

J ; = L(log Pt-*/) + r(log Pt-*/). (6.52) 

-TT-tJ 

For such /'s we now introduce the two functionals 

<I>i(t) =Pa(Pr-J)r(lnPT-t/)), 

<I>2(t) = Pt {{Pr-t f )r'' {In Pr-tf)) . 

Notice that for every T > the function u{x,t) = Px-tfix) satisfies the backward Cauchy 
problem 

(Lu + ^ = u > in M X (-oo,r). 



dt 

u{x,T) = f{x). 



It follows that 



|<I.i(0) = PT(/)r(lnPT/), <1>2(0) = PT(/)r^ (In Pt/), 
\$i(r) = PT(/r(ln/)), ^2{T) = PrifT^ilnf)). 

Lemma 6.4 We have 

<!>[{t) = 2Pt ((Pr_t/)r2(lnPT^ J)) , 

and 

$'2(t) = 2Pt ((Pr-J)rf (InPr-i/)) . 



(6.53) 
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Proof. Let ^ be a smooth vector field on M and let $y be the functional 

$y(i) = Pt {{PT-tf){VlnPT^tff) . 

We have 

<!>'y{t) = Pt {L{{PT-tf){V In PT-tff)) - Pt {{LPT-tf){V In Pr-tff) 

'LPr^tf 



We now compute 



2Pt[{VPT-tf)V , 



Pt {Li{PT-tf){V In Pr-tff)) = Pt {L{PT-tf){V\nPT^tff)+Pt {{PT^tf)L{V\nPT-tff) 
+ 2Pt{T{PT-tf, {V\n PT-tff)) 

= Pt {L{PT-tf){V In Pr-tf?) + 2Pt {{Pt -t f ) {V In Pt -tf){LV In Pr-tf)) 

+ 2Pt {{PT-t f )r{V In PT-tf,V In PT-tf))+4Pt {V In PT-tfT {PT-tf,V In PT-tf))) ■ 

By taking into account (j6.52p we obtain 

<^>'y{t) =2Pt {{V PT-tf )[L, V] In PT-tf) + 2Pt {{PT-tf)r{V In PT-tf,V In PT-tf)) 

+ iPt {V In PT-tfT{PT-tf,V In PT-tf))) - 2Pt {{V PT-t f )Vr {In PT-tf, In PT-tf)) ■ 

We now observe that 

d 

Vr{ln PT-tf, In PT-tf) = 2T{ln PT-tf, Z In PT-tf) + 2Y,{Xi\n PT-tf ){[Z,Xi] In PT-tf). 

1=1 

Thus, 

<!>'y{t) = 2Pt {{PT-tf)rX {In PT-tf, In PT-tf)) (6.54) 
-APt ^{VPT-tf)J2{XilnPT-tf){[V,Xi]lnPT-tf)^ , 

where we have defined 

r^(/,5) = \{L{{Vf){Vg)) - VfVLg - VgVLf. 
We first apply (I6.54p with V = Xj, and sum in j = 1, d to obtain 

<I>;(t) = 2Pt {{PT-t f )r2{ln PT-tf)) . 
Here, we have used to skew-symmetry of which gives 

d d 

Y,{XjPT-tf)Y,{X,ln PT-tf ){[X,,Xi]ln PT-tf) = 0. 
j=i i=i 



37 



We next apply (j6.54p with V = Zmn and sum in m, n = 1, f), obtaining 
$'2(t) = 2Pt {{Pr-t f )rl {In Pr-tf)) 

-4 ^ Pt({Z^nPT-tf)Y.{X,lnPT-tf){[Z^n,X,]lnPT-tf)\ 
m,n=l \ i=l / 

We now observe that, thanks to the crucial assumption ()2.5p . one has 
f) d 

J2 {ZmnPr-tf) Y^^Xi In PT^tf){[Zran, X,] In Pr-tf) 
m,n=l 1=1 

t) d 

= - Yl {ZmnPT-tf) Yl SLniXiln PT-tf){XelnPT-tf) = 0. 

m,n=l i,£=l 

Therefore, 

$'2(0 = 2Pt ((PT-J)rf (lnPT_i/)) . 



□ 



Proposition 6.5 Let b be a smooth, positive and decreasing function on the time interval [0, T]. 
On M X [0,r], we have 

-— d>i + 6$2 >--r-^LPTf + -^PTf, 
2/02 / ap2 dp2 

where 

d fb" Kb' ^ \ 

Proof. To prove this result we apply the sub-Riemannian curvature-dimension inequality (|6.46|) , 
in combination with Lemma |6.4[ If a and b are positive functions we thus obtain 

(a$i + 5$2)' > [a + 2pia - 2k^^ «>i + {b' + 2p2a)^2 + ^ (Pt((PT- J)(ilnPT-t/)')) 
But, for every 7 G R, 

(LlnPT_J)2 >27LlnPT-J-7', 

and 



L In Pr-tf = - r(ln Pr-tf). 



Therefore, 



(a$i + &$2)' > fa' + 2pia - 2k^ - 4^^ $1 + (6' + 2/92a)$2 + ^LPrf - 2^Pt/. 

\ b a J a a 

By taking a, 6, 7 such that 

a2 4 

a + 2pia — 2k- -70 = 0, 

b d 

b' + 2p2a = 0, 

we obtain the desired conclusion. □ 
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6.3 Li-Yau type estimates 

In this section, we extend the celebrated Li-Yau inequahty in [IB] to the heat semigroup asso- 
ciated with the subehiptic operator L. Let us mention that, in this setting, related inequalities 
were obtained by Cao-Yau [13] • However, these authors work only in the case of a compact 
manifold and do not base their study on the analysis of a tensor like TZ. As a consequence they 
do not obtain a control of the constants in terms of the geometry of the manifold. 

Proposition 6.6 (Gradient estimate) Assume that (|6.45p hold. Let f e Co°(M), with f > 
0, then the following inequality holds for t > 0; 



r(inPj) + -^tr^(inPj)< 1 + - ^ + 1+ + 



2 



2/02 ^ J Ptf 6 2 V 2p2j 2t 
Proof. If we apply Proposition 16.51 with b{t) = (T — t)^, we obtain: 

d f2piiT-t)-^-2\ d f 1 A , 3k 



and 



The result easily follows. □ 

Remark 6.7 When pi > 0, then (j6.45p and (j6.46p a/so /loW in particular with pi = 0. a 
consequence of Proposition 16.61 



r(lnPJ) + ^tr^(lnPJ) < (1 + 1^)^ + lii^i^. (6.55) 
3 2p2 Ptf 2t 

This inequality leads to a optimal Harnack inequality only when pi = (which corresponds to 
Ric> 0). Sharper bounds in the strictly positive curvature case will be obtained in (112. 89p by a 
different choice of the function b{t). 

Remark 6.8 // we set 

« = d(l + |^). (6.56) 

then, as a consequence of (|6.55p . we obtain that in the case pi > 0, 

-Pj - -2t ^'-''^ 

The constant — y m (j6.57p is, in general, not sharp, as the example of the heat semigroup 
on a graded nilpotent Lie group shows. In such case, in fact, one can argue as in l2W to 
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show that the heat kernel p{x,y,t) is homogeneous of degree —-j with respect to the parabolic 
dilations {x,t) {5{X){x), X'^t), where (5(A) represent the non-isotropic dilations associated with 
the grading of the Lie algebra of M, and Q = d + 2\] indicates the corresponding homogeneous 
dimension ofM.. From such homogeneity ofp{x,y,t), a scaling argument produces the estimate 

LPtf > _Q 
Ptf - 2t' 

which, unlike (|6.57|) . is best possible. However, the estimate (|6.57|) is sharp in the case of the 
Laplace operator in for which we have k = 0. It seems difficult to obtain sharp geometric 
constants by using only the inequality (j6.46p . In part this is due to the fact that when we apply 
Cauchy-Schwarz inequality in the Bochner formulas, we loose a piece of information. Another 
difference is in the presence of the commutator terms in (]6.46p . These aspects are quite different 
from the Riemannian case, for which the CD{d,R) inequality (|l.ip provides sharp geometric 
constants (see J^, 

7 Stochastic completeness of the heat semigroup 

In [75] Yau proved that if M is a complete Riemannian manifold with a Ric > p, with /? G R, 
then one has the stochastic completeness of the heat semigroup, i.e. Pfl = 1. Under the same 
assumptions, Dodziuk [23] proved that bounded solutions of the heat equation are characterized 
by their initial condition. 

In this section, we extend these theorems of Yau and Dodziuk to our setting. Throughout this 
section, we assume that (M, d) is complete and that for every smooth function / : M ^ M and 
every > 0, 

T2{f, f) + uTUf, f) > \{Lff + (pi - ^) r(/, /) + p2r^(/, /), (7.58) 

where pi G M, whereas p2,n > As in the previous section, the idea is to study monotone 
increasing functionals of the heat semigroup. 

Proposition 7.1 Let f G C^(M) and T > 0. If a,b : [0,T] M are two positive functions 
such that 

a'{t) + 2pia{t) - > 

b'it) + 2p2a{t) > 0, 

then the functional 

E{t) = a{t)Pt{T{PT.tf)) + h{t)Pt{T^{PT-tf)) 

is monotone increasing. 

Proof. We begin by noting that, with <^i(t),$2(i) as in section W72[ we have 

E{t) = a{t)^i{t) +b{t)^2{t). 
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Lemma 16.41 now gives 

E'it) = a'{t)Pt{T{PT-tf)) + h'{t)Pt{T^{PT-tf)) + 2a{t)Pt{T2{PT^tf)) + 2b{t)Pt{Tl{PT-tf)). 
Using the curvature-dimension inequality (j6.46p we now find 

a{t)Pt{T2{PT^tf)) + b{t)Pt{Ti{PT^tf)) 



Therefore, 

E'{t) > (^a'(t) + 2pMt) - ^^^) PtinPT-tf)) + {b'{t) + 2p2a{t))Pt{T^{PT^tf)) > 0. 
The desired conclusion immediately follows from this inequality. 

□ 

The following energy estimate represents a basic consequence of Proposition 17.11 

Corollary 7.2 There exists v ^ (v < 2min{/>2,pi — will do), such that for every f G 
C^(M), one has 

r{Ptf) + T^iPf) < e-^* (PtTif) + PtF^if)) . (7.59) 
As a consequence, for every 1 < p < oo one obtains 

l|r(Pt/)llLP(M) <e-'^*(l|r(/)||LP(M) + l|r^(/)llLP{M)), t>o. (7.60) 

Proof. Let t > and consider the interval < s < t. With the choice 

a(s) = b{s) = e'"', 0<s<t, 

where z/ € M satisfies v < 2min{/32,pi — k}, Proposition 17.11 gives £'(0) < E{s). This inequality 
reads 

r(pj) + r^(pj) < e-^' {PsiriPt-sf)) + Psir^iPt^sf))) • 

Letting s ^ t~ we obtain (j7.59p . Combining (|7.59p with (16.49P we deduce (|7.60p . 

□ 

Theorem 7.3 For t > 0, one has Ptl = 1. 

Proof Let f,g£ C^(M), we have 

/ {Ptf - f)gdfi = f [ (-^Psf)gdfids= f I iLPsf)gdfids = - f [ r{Psf,g)dfids. 
Jm Jo Jm J Jo Jm Jo Jm 

By means of ()7.60p . and Cauchy-Schwarz inequality, we find 

/ {Ptf-f)gdp < ( fe-'^ds) J||r(/)||oo + ||r^(/)|U / r{g)hdp. (7.61) 

Jm \Jo / JM 
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We now apply (j7.6ip with / = where are the functions in Lemma 16.11 and then let 
n — > oo. Since by Beppo Levi's monotone convergence theorem we have Pthn{x) /" Ptl(x) for 
every x G M, we see that the left-hand side converges to jy^{Ptl — l)gdfj.. We claim that the 
right-hand side converges to zero. To see this observe that, thanks to the assumption (j2.1ip . we 
have for any cj) G C°°(M) 

Therefore, 

^^(MU^S^^(MU <C||Vfi/i„||oo ^0, asn^cx). 
We thus reach the conclusion 

/ {Ptl - l)gdfi = 0, geC^{M). 
Jm 

It follows that Ptl = 1- 

□ 

Theorem 17.31 is equivalent to the uniqueness in the Cauchy problem for initial data in L°^(M). 
Proposition 7.4 The unique bounded solution of the Cauchy problem 

^-Lu = 0, 

uix,0) = fix), /GL-(M), 
is given by u{x,t) = Ptf{x). 

Proof. Since the vector fields Xj's are locally Lipschitz (because smooth), for every x G M, the 
stochastic differential equation 

= x+ f Xo{Y^)ds + V2y] fxi{Yf)odBi, 
Jo Jo 

has a solution up to a stopping time e(x), where {Bt)t>o is a d-dimensional Brownian motion. 
It is seen that 

{Qtf){x) = E(li<e(,)/(11")), / G L\M), 

defines a contraction semigroup on L^(M). By Ito's formula its generator on Co(M) is given by 
L. By uniqueness of the heat semigroup, we actually have Qtf{x) = Ptf{x), for / G Co(M). By 
using the definition of Qt, we deduce that for / G Co(M), 

Pt/(x)=E(lt<e(.)/(yt")) 

By applying the previous equality with fhn, and letting n — > -|-oo, we deduce from Beppo Levi's 
monotone convergence theorem that 

Pt/(x) = E(ii<e(,)/(y,-)) 
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when / is a bounded and smooth function. Since Pfl = 1, almost surely e{x) = +00, and 

Ptf{x) = E{f{Yf)). 
Let now u{t, x) be a bounded solution of the Cauchy problem: 

\u{x,0) = fix), 

where / is a bounded and smooth function. From Ito's formula, for t > the process {u{t — 
s, Yf))o<s<t is a local martingale and since u is bounded, it is a martingale. The expectation of 
{u{t — s, Y'/))o<s<t is therefore constant. From this we deduce 

u{t,x)=Eif{Yf)) = Ptf{x). 

Finally, when / E L°°(M), for the Cauchy problem 



du 



Lu = 0, 



[u{x,0)=fix) 

we have from the regularization property of the heat semigroup and the above uniqueness result: 
For every r > 0, 

u{x,t + T)= / p{x,y,t)u{y,T)fj.{dy). 



By letting r — > 0, we obtain u{x,t) = Ptf{x) as desired. □ 



8 A parabolic Harnack inequality 

In this section we generalize the celebrated Harnack inequality in [36] to nonnegative solutions 
of the heat equation H = L — on M. which are in the form u{x,t) = Ptf{x), for some 
/ E C°°(M) n L°°(M). Theorem 18.11 below should be seen as a generalization of (i) of Theorem 
2.2 in [46], in the case of a zero potential q. One should also see the paper [14J, where the 
authors deal with subelliptic operators on a compact manifold. As we have mentioned, they do 
not obtain bounds which explicitly depend solely on the geometry of the underlying manifold. 

Theorem 8.1 Let M. be a complete sub-Riemannian manifold such that (|6.45|) holds with pi > 
0. Let f E C°°(M) be such that < f < M, and consider u{x,t) = Ptf{x). For every 
(x, s), (y, t) E M X (0, 00) with s < t one has with D as in (16.560 

„(...,<„fa,*)0)*exp(f||^). (8.02) 

Proof. Let / be as in the statement of the theorem, and consider /„ = hnf, where /i„ is the 
sequence in Lemma [6.11 Then, fn E Cq°{M), and < fn y /• By Beppo Levi's monotone 
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convergence theorem we have n„(x, t) = Ptfn{x) /" ii(x, t) = Ptf{x) for every (x, t) £ Mx (0, oo). 

at 



Since Lm„ = ^§r, in terms of ii„ inequality (j6.55p can be reformulated as 



A / /I , 3k 9 log n„ ^(^+2^2) 



r(in n„) + -^tr^ (In + , 



+ 



2i 



In particular, this implies 



2t 



(8.63) 



We now fix two points (x, s), {y, t) S M x (0, 00), with s < t. Let 7(t), < r < T be a subunit 
path such that 7(0) = y, 7(T) = x. Consider the path in M x (0, 00) defined by 

a(r) = (7(r),t + ^r ) , < r < T, 



so that a(0) = {y,t), a(T) = {x,s). We have 



In 



Unix,s) _ r _d_ 

Uniy,t) Jo (It 

■T 



In Un{a{T))dT 



f — s (/ In 7/ 

< 7'(r), V^Xlnn„)(a(r)) > -— ^(a(T)) 



T 5t 

Now since 7(t) is subunitary we have 

|< 7'(r), V/j(lnu„)(a(r)) >| < r(lnu„(a(r)))^ 
and applying ()8.63p for any e > we find 



dT 



log 



Unjx.s) 

Un{y,t) 



'0 

< — r + - 

- 2e 2 



2P2 7 



If we now choose e > such that 



r(ln Un){a{T))dT 

T 

r(ln Un){a{T))dT 



dT 



t - s 



(91n ; 



dt 



— {a{T))dT 



t - s 



r(ln Un){a{T))dT 



t+ -j^T 



e t — s 

2 = r(i + ^)' 



we obtain from the latter inequality 

..(x,.)^ 4(7)^(1 + ^) , + 



log 



Un{y,t) 4(f-s) 



+ 



In 
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where we have denoted by ^^(7) the subunitary length of 7. If we now minimize over all 
subunitary paths joining y to x, and we exponentiate, we obtain 

, ,^ , ./aK^+H) (d{x,y)Hl + $-j\ 
Un{x,s) <Un{y,t) [ - ] exp ' 



Ait - s) 



Letting n — > 00 in this inequality we finally obtain (j8.62p . 



□ 



The following result represents an important consequence of Theorem 18.1 



Corollary 8.2 Let p{x, y, t) be the heat kernel on M. For every x,y,z£Wl and every < s < 
t < 00 one has 



p{x,y,s) < p{x,z,t) - exp 



t\^ (Dd{y,zf 



sj \d A{t- s) 

Proof. Let r > and x G M be fixed. By Hormander's hypoellipticity theorem [37] we know 
that p{x, •, • + r) G C°°(M X (-r, 00)). From (f630|) we have 

p{x, y,s + T) = Ps{p{x, ■,T)){y) 

and 

p{x,Z,t + T) = Pt{p{x,-,T)){z) 

Since we cannot apply Theorem 18.11 directly to u{y,t) = Pt{p{x,-,T)){y), we consider as in the 
proof of the latter Un{y,t) = Pt{hnp{x, •,r))(y), where hn G Cg°(M), < /i„ < 1, and /i„ 1- 
Prom ([8:62]) we find 

t\^ (Dd{y,z) 



Ps{hnp{x,-,T)){y) < Pt{hnp{x,-,T)){z) - exp 



d A{t - s] 

Letting n — > cxd, by Beppo Levi's monotone convergence theorem we obtain 



p{x, y,s + t) < p{x, z,t + t) [ - \ exp 



t\^ (Dd{y,zf 



The desired conclusion follows by letting r — > 0. 



n 



9 OfF-diagonal Gaussian upper bounds for p{x,y,t) 

Let M be a complete sub-Riemannian manifold such that (j6.45p holds with pi > 0. Fix x G M 
and t > 0. Applying Corollarv 18.21 to {y,t) p{x,y,t) for every y G B{x,\/t) we find 

p{x, x, t) < 2^635 p{x, y, 2t) = C{d, k, P2)p{x, y, 2t). 
Integration over B{x, \/t) gives 

p{x, X, t)p.{B{x, Vt)) < C{d, K, P2) / p{x, y, 2t)dp{y) < C{d, k, P2), 
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where we have used Ptl < 1. This gives the on-diagonal upper bound 

The aim of this section is to estabhsh the following off-diagonal upper bound for the heat kernel. 

Theorem 9.1 Let M. be a complete sub-Riemannian manifold such that ()6.45p holds with pi > 
0. For any < e < 1 there exists a constant C{d, k, p2,€) > 0, which tends to oo as e — > 0"*", 
such that for every x, y € M and t > one has 

p{x,y,t)< C{d,K,p„e) ^^^^f d{x,yf 



p{Bix,Vi))-^p{B{y,Vt))-^ V (4 + e)t 

Proof. We suitably adapt here an idea in [14] for the case of a compact manifold without 
boundary. Since, however, we allow the manifold M to be non-compact, we need to take care 
of this aspect. Corollarv 17.21 will prove crucial in this connection. Given T > 0, and a > we 
fix < r < (1 a)T. For a function ip £ C^(M), with ijj > 0, inMx (0,r) we consider the 
function 

fiy,t)= Piy,z,t)p{x,z,T)'ilj{z)dp{z), X G M. 



Since / = Pt{p{x, •,r)V'), it satisfies the Cauchy problem 

iLf-ft = inMx(0,r), 
\fiz,0)=p{x,z,T)ij{z), zGM. 

Notice that thanks to Hormander's theorem [37] we know y — > p{x,y,T) is in C°°(M), and 
therefore •,r)V' G L°°(M). Moreover, (j6.49p gives 



||Pt(p(x,-,T)V')||i2(M) < ||p(x,-,r)V'||i2(M) ^ J^P^^'^''^^ i^{z)dp{z) < oo, 

and therefore ^ 

/ / f{y,tfdp{z)dt <T p{x, z,Tf7p{z)dp{z)dt < oo. (9.65) 
Jo Jm Jm 

Invoking ()7.59p in Corollary 17.21 we have 

r{f){z,t) < e-"' {PtT{p{x,;T)^){z)+PtT^{p{x,;T)^){z)) . 
This allows to conclude 

T{f){z,tfdp{z)dt<e\^\- [ {r (pix, ■,T)^)izf +T^{pix,;T)^)izf } dp{z) < oo. 

Jm 

(9.66) 

We now consider a function g G C^([0, (1 + a)r], Lip^(M)) n L°°(M x (0, (1 + a)T)) such that 

-^>^r(5), onMx (0,(1 + Q)r). (9.67) 
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Since 



(L - = 2/(L - |)/ + 2r(/) = 2r(/), 



multiplying this identity by h^{y)e^^'^'^\ where hn is the sequence in Lemma [6. H and integrating 
by parts, we obtain 



= 2 
= 2 
+ 2 



hle3T{f )d^,{y)dt - 
hle<^T{f)d^Ji{y)dt + 4 



hle^fT{f,g)dfxiy)dt- f [ hnS^ f%f,{y)dt 



hle^{L-^)fdpi{y)dt 
hne^ fT{hn, f)dn{y)dt 



Ke>^fdfi{y) 



+ / Ke^fdiiiy) 



t=T 



>2 r [ hleS V fx,/ + l-Xig\ dfi{y)dt + 4 T [ Ke^ fT{K, 
Jo Jm i=i V ^ / Jo Jm 

+ / Kesfdfiiy) - [ KeSfdfiiy) 
Jm t=T Jm 



f)dii{y)dt 



where in the last inequality we have made use of the assumption (j9.67p on g. From this we 
conclude 



Ke3fd^i{y) 
We now claim that 



t=T 



< / Ke^^fdfiiy) 



t=o JO 



Ke^ fT{hnJ)dii{y)dt. 



lim 

n^oo 



hne3fT{KJ)dii{y)dt = ^. 
To see this we apply Cauchy-Schwarz inequality which gives 



< 



Ke^ fT{hn, f)dfi{y)dt 
eafT{hn)dfi{y)dt 



< 



hle^fr{hn)dfi{y)dt] / / ear{f)dfi{y)dt 



e3T{f)dfi{y)dt\ ^0 



as n — > oo, thanks to (j9.65p . (j9.66p . With the claim in hands we now let n — > oo in the above 
inequality obtaining 



eSiy<-)f\y,T)df,iy)< / e^(^'0)/2(y,0)(i/i(y). 



(9.68) 



At this point we fix j; G M and for < t < r consider the indicator function of the ball 

B{x,y/t). Let ipk G C'o°(^)' i'k > 0, be a sequence such that ipk '^B{xy/i) L'^{M.), with 
supp ipk C B{x, 100\/i). Slightly abusing the notation we now set 



fiy,s) = Ps{p{x,-,T)lj^^^^^^){y)= / p{y,z,s)p{x,z,T)dfi{z). 
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Thanks to the symmetry of p{x, y, s) = p{y, x, s), we have 



p{x,z,Tfdfi{z). 



'B(x,Vt) 

Applying (|9.68p to fk{y,s) = Ps{p{x, ■,T)^l)k){y), we find 

Jm 

At this point we observe that as /c — > oo 

JM 

< 2||e^'(''^)||Loo(M)||p(x,-,r)||i2(M)||p(x,-,r)||^oo(B(^^iio7t))||^fc - 1^(^,^4)1 1 l2(m) 



(9.69) 



(9.70) 



0. 



By similar considerations we find 

3^(^'°Vl(y,0)dM(y) 



5s(^'0)/2(y,0)dM(y) 



< 2||e3(''°)||2.oo(M)||p(2;,-,r)||^^(^(^^^^o^))||V'fc - lB{x,Vi)llL2(M) 



0. 



Letting /c ^ cxd in (j9.70p we thus conclude that the same inequality holds with fk replaced by 
f{y, s) = Ps{p{x, -jT)!^^^ This implies in particular the basic estimate 



inf e»(^'^) 

z(^B{x,Vt) 



f\z,T)dll{z) 



< 



B(x,y/t) 

3^''-^f{z,T)df,{z)< I e^^''^^f{z,T)df,{z) 



(9.71) 



IB{x,Vt) 

< [ e^(^'0)/2(z,0)d^(z) 



< sup e 

zeB(y,Vi) 

At this point we choose in ()9.7ip 



9iy,t) = 9x{y,t) 



'B(y,Vt) 

p{x,z,Tfdfi{z). 

d{x,yf 



?^''''^^p{x,z,Tfdfi{z) 



2((1 + 2a)r-t) 



Using the fact that T{d) < 1, one can easily check that (j9.67p is satisfied for this g. Taking into 
account that 

inf (^9x(Z,T) _ g 2{{1+2c«)T-t) > g2{{l+2a)T-T) ^ 

zGB{x,Vt) zeB{x,Vi) 

if we now choose r = (1 + a)T, then from the previous inequality and from (|9.69p we conclude 
that 



f\z, (1 + a)T)dfi{z) < sup e WT^' 

B{x,Vt) \zeB{y,Vt) 



B(y,Vt) 



p{x,z,Tfdfi{z). (9.72) 
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We now apply Theorem 18.11 which gives for every z G -B(x, 



Integrating this inequahty on -B(x, \ft) we find 



/ p{x,z,Tfdti{z) 



3« N 



,2 {\^af^^^^\^^ I „2, 



fix, TY < ^ / f\z, {l+a)T)dfi{z). 

IJ.[B{x,Vt)) JB{x,Vi) 



If we now use (|9.72p in the last inequality we obtain 

/ p{x,z,Tfd^x{z) < ^ ^— sup e 2(i+2„)T^2.T . 

JB{y,Vt) fl{B{x,Vt)) \zeB{y,Vt) J 

Choosing T = (1 + a)t in this inequality we find 

^ (1+ •) ^ 
(1 _|_ Q,')'^(^+2^)g2c«(l + a) +2a(l + Q) / d(a;,z)2 1 

(1 +Q)t)^d/^(Z) < ^ — sup e 2{l + 2a)(l + <.)t + 2a(l + a) 

\zeBiy,Vi) 

(9.73) 



We now apply Corollary 18.21 obtaining for every z £ B{y, ^/t) 



1 + — 

^ ^ 2p2 

2a 



p{x, y, tf < p{x, z, (1 + a)tf{l + a) V +2p2; exp 
Integrating this inequality in z B{y, \/t), we have 

n{B{y, Vi))p{x, y, tf < (1 + af^^^^) [ p{x, z, (1 + a)tfdn{z). 



B{y,Vt) 



Combining this inequality with (j9.73p we conclude 



p{x,y,t)<- J 1 sup e 2(i+2<.)(i+c.)t 

fi{B{x,Vt))2n{B{y,Vi))2 \z(LB(y,Vt) 

If now x G B{y, \/t), then 

z)2 > {d{x, y) - Vtf > d{x, yf - t, 

and therefore 



sup e 2{l + 2a)(l + a)t < g 2(1 + 2q) + g 2(l + 2a){l + a)f _ 

If instead x ^ B{y, then for every 5 > we have 

d{x, zf > (1 - 6)d{x, y f - (1 + (5~^)t 
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Choosing (5 = q/(q + 1) we find 



d{x,zf > 



{2 + a-^)t, 



l + a 



and therefore 



g 2(l + 2a){l + c«)t < g 2{l+2a)(l + a)^t 2(l + 2a){l + a) 




sup 

z&B(y,Vi) 



For any e > we now choose a > such that 2(1 + 2a) (1 + a)^ = 4 + e to reach the desired 



10 A generalization of Yau's Liouville theorem 

In his seminal 1975 paper [73], by using gradient estimates, Yau proved his celebrated Liouville 
theorem that there exists no non-constant positive harmonic function on a complete Riemannian 
manifold with non-negative Ricci curvature. The aim of this section is to extend Yau's theorem 
to the sub-Riemannian setting of this paper. An interesting point to keep in mind here is 
that, even in the Riemannian setting, our approach gives a new proof of Yau's theorem which 
is not based on delicate tools from Riemann geometry such as the Hessian and the Laplacian 
comparison theorems for the geodesic distance. However, due to the nature of our proof at the 
moment we are only able to deal with harmonic functions bounded from two sides, whereas in 
|73j the author is able to treat functions satisfying a one-side bound. 

In what follows we assume that the curvature-dimension inequality (j6.46p hold with pi = 0, 
P2 > and «; > 0. We also assume that the metric space (M, d), or equivalently (M, d/j), is 
complete. As we mentioned it before, see Remark 15.31 the Riemannian case of Example 12.11 
our assumptions are in fact equivalent to assuming that M is a complete Riemannian manifold 
with Ric> 0. We begin with a Harnack type inequality for the operator L. 

Theorem 10.1 Let M 6e a complete sub-Riemannian manifold and assume ()6.46p with pi = 0. 
Let < f < M be a harmonic function on M, then there exists a constant C = C{d, p2, > 
such that for any xq gM. and any r > one has 



Proof. Let / be as in the statement of the theorem. By Hormander's theorem j37j we know 
that / G C°°(M). Applying Theorem 18.11 to the function u{x,t) = Ptf{x), we obtain for 

x,y £ B{xo,r) 



At this point we observe that, thanks to the assumption Lf = 0, the functions u{x,t) = Ptf{x) 
and v{x,t) = f{x) solve the same Cauchy problem on M. By Proposition 17.41 we must have 



conclusion. 



□ 



sup f < C inf /. 

B{xo,r) B(xo,r) 




< s < t < oc. 
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Ptf{x) = f{x) for every x G M and every t > 0. Therefore, taking s = r'^,t = 2r^, the latter 
inequahty gives 

/(x) < (\/2es) /(y), x,y e B{xo,r). 

n 

Corollary 10.2 (of Cauchy-Liouville type) Assume (j6.46p with pi = 0, then there exist no 
bounded solutions to Lf = on M, other than the constants. 

Proof. Suppose a < / < 6 on M. Consider the function g = /— inf /. Clearly, < g < M = b—a. 

M 

If we apply Theorem 110.11 to g we find for any xq G M and r > 

sup g < C inf g. 

B{xo,r) B{xo,r) 



Letting r — > oo we reach the conclusion sup / = inf /, hence / = const. 

M M 



□ 



11 Volume growth and Isoperimetric inequality when pi = 

Throughout this section, we shall assume that (M, d) is complete non compact and that (|6.45p 
holds with pi > 0. 

11.1 Volume growth 

We first derive a basic and straightforward consequence of the parabolic Harnack inequality on 
the volume growth of metric balls. 

Proposition 11.1 For every x G M and every Rq > there is a constant C{d, k, P2) > such 
that, with D as in ()6.56p . 



/x(B(x,i?))<-g^f^l^ii^, R>Ro. 
R^p{x,x,R^) 

Proof. Let t > r > 0. From the Harnack inequality of Corollarv 18.21 we have 

n 

p{x,x,t) > p{x,x,t) i^- 



On the other hand, the inequality (j9.64p gives 

C{d,K,p2) 



p{x, X, t) < 



p{B{x,Vi)) 
This implies the desired conclusion. 



□ 
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11.2 Isoperimetric inequality 

In [31] it was proved that in a Carnot-Caratheodory space [X, ^, d) the doubling condition 

^^{B{x, 2r)) < Cifi{B{x, r)), x e X,r > 0, 

for the volume of the metric balls combined with a weak Poincare inequality suffice to establish 
the following basic relative isoperimetric inequality 



mmMEnB{x,r)),fx{{X\E)nB{x,r))}— <aso[^, ^ P{E,B{x,r)), (11.74) 

where P{E, B{x,r)) represents a generalization of De Giorgi's variational notion of perimeter, 
and E C X is any set of locally finite perimeter. In this inequality the number D = log2Ci, 
where Ci is the doubling constant, and Cjso is a constant which depends only on Ci and on the 
constant in the Poincare inequality. If in addition the space X satisfies the maximum volume 
growth 

l^{B{x, r)) > C2r^, x G M, r > 0, (11.75) 

then (|11.74p gives 

mm {n{E n B{x,r)), ^i{{X \ E) n B{x,r))}^ < C*,^P{E, B{x,r)), (11.76) 
where C*„ = C^'^^'aso- 

When X = M is a Riemannian manifold with Ric > 0, then the doubling condition with 
Ci = 2^ , D = dim M, follows from Bishop-Gromov comparison theorem (see [19], Proposition 
3.3 and Theorem 3.10), whereas the Poincare inequality was proved by Buser, see [13J. As a 
consequence, one obtains the relative isoperimetric inequality (jll.74p in this setting. When M 
satisfies the maximum volume growth (jll.75p . one also obtains from (|11.74p the global isoperi- 
metric inequality 

fi{E)^ <C\,oP{E,m), (11-77) 
for any measurable set of locally finite perimeter E C M. 

In this subsection we investigate the sub-Riemannian counterpart of the isoperimetric estimate 
()11.77p under the assumption that pi = 0. Here, the main obstacle is precisely the a priori lack 
of a global doubling condition and of a Poincare inequality, and therefore we cannot rely on the 
above cited results from [31]. Instead, using our Li-Yau type estimate (I11.79P we adapt some 
beautiful ideas of Varopoulos and Ledoux to provide a characterization of those sub-Riemannian 
manifolds which support an inequality such as (jll.77p . We stress that, due to the nature of our 
approach, we obtain a lower bound on the dimension D in (lll.77p . but in the case of a graded 
nilpotent Lie group our D is not optimal, see Remark 16.81 
In what follows, given an open set C M we will indicate with 

= {(p e c^{n,n) \ ||</.||oo < i}. 
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Here, for cp = Yli=i 4>i^i^ we have let llc/illoo = sup 

yEi=i<Pl Following Ll5j, given a function 
n ^ 

f G Lj^^{i^) we define the horizontal total variation of / in as 




The space 

BVnm = {/ e L^n) I Var7^(/; Q) < cx,}, 

endowed with the norm 

\\f\\BVn{n) = II/IIli(m) + Var7^(/;n), 

is a Banach space. It is well-known that W^^{Q) = {/ € L^{Q,) \ Xif G L^{^),i = l,...,d} is a 
strict subspace of BV'^{Q,). It is important to note that when / G W^^{Q), then / € BV-)^{^), 
and one has in fact 

Var^(/;17) = ||7f(7)||Li(n). 

Given a measurable set S C M we say that it has finite horizontal perimeter in 17 if 1^; G V?^ (17) . 
In such case the horizontal perimeter of E relative to is by definition 

We say that a measurable set C M is a Caccioppoli set if Pn{E; Q) < oo for any d M. We 
will need the following approximation result, see Theorem 1.14 in [31j. 

Proposition 11.2 Let f € BV-ni^), then there exists a sequence {/njneN of functions in 
C°°{n) such that: 

(^) ll/n-/||Li(f7)-0; 

(ii) VT{f;:)dfi ^ Varnif; n). 

If Q = M, then the sequence {fn}neN ca7^ be taken in C^(M). 

Our intent is to establish the following result. 

Theorem 11.3 (Isoperimetric inequality) Suppose that M is not compact in the metric 
topology and that there exists D > 1 such that 

fi{B{x,r))>Cir^. (11.78) 

Then, there is a constant Cjso = Ciso{d, P2, i^,Ci, D) > 0, such that for every Caccioppoli set 
E CM 

KE)^ <CisoPn{E,M)- 

The essence of the proof of Theorem 111.31 is contained in the following result. 

Proposition 11.4 Let D > 1. Let us assume that M is not compact in the metric topology, 
then the following assertions are equivalent: 
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(1) There exists a constant Ci > such that for every x £ M, r > 0, 

n{B{x,r)) > Cir^. 

(2) There exists a constant C2 > such that for x G M, t > 0, 

p{x,x,t) < 

t^ 

(3) There exists a constant C3 > such that for every Caccioppoli set E cWl one has 

fi{E)^ <C3Pn{E;M)- 

(4) With the same constant C3 > as in (3), for every f G BVni^) ^^-^ 

D-l 

Proof. That (1) (2) follows immediately from the Gaussian bound of Theorem 19.11 (observe 
that we may take C2 = ^). 

The proof that (2) (3) is not straightforward. First we note that (j6.45p implies the Li-Yau 
type estimate (I6.55p . This enables us to adapt some beautiful ideas of Varopoulos (see [70] . 
pp.256-58) and Ledoux (see pp. 22 in [43j, see also Theorem 8.4 in [44])- Let / G Co(M) with 
/ > 0. By ([635]) we obtain 

mf) - (1 + l^w^ < ^-^^{Ptff. (11.79) 

This gives in particular, with . = ^ (l + ^ 

^"^'^^^ < ^Ptf, (11.80) 



dt J - 2t 

where we have denoted = sup{a, 0}, a~ = sup{— a,0}. It follows that for every < Ti < 
T2 < 00 

V L{^) ^/^^*^il"(^2/Ti)||/||ii(M)<oo. (11.81) 
Using Tonelli's theorem we now have for any t > 



Ptf{x)dfi{x) = / f{y) / p{x,y,t)dfi{x)dfi{y) = / f{y)dfi{y), 
; Jm Jm Jm 

where we have used Theorem 17.31 This gives for any t,h > 

= [ Pt+hf{x)dfiix)- [ Ptfix)dfiix) = [ f^^ ^^(x)dsd^iix) 
Jm Jm Jm Jt os 

—^{x)dfi{x)ds, 



54 



where the exchange of order of integration is justified by Fubini's theorem, which we can apply 
in view of (jll.Sip . The latter equation implies 



t+h 



dfi{x)ds = 



dPsf, 



and therefore, we have 



t+h 



dPsf 



ds 



x) 



d^{x)ds 



< 



t+h 



(x) I dfi{x)ds, t,h>0, 



(x) I dfi{x)ds 



dPsf, 
ds 



ds_ 

s 



Li(M)- 



Letting /i ^ 0+ we finally conclude 



I|-^IIli(M) ^ yll/llLi(M)' * > 0. 

By duality, we deduce that for every / G Co(M), / > 0, 



^dPtf, 



Il°°(m) 



< 



Once we have this crucial information we can return to (jll.79p and infer 



r(Pt/)<^^ll/llioo(M), i>o. 



Thus, 



ll\/r(Pt/)||ioo(M) < -^\\J ||L°°(M)- 

Applying this inequality to 5 € C^(M), with g > and \\g\\L°°{M) ^ f ^ C'o(M) we have 

/ 9{f - Ptf)dn = [ [ g^^dfids = 
Jm Jo Jm OS 



PsLgfd^ds 



< 



wVnPs 



9)\\l^ 



gLPsfd^ds = LgPsfd^ds 
Jo Jm 

LPsgfdfids = - / T{Psg,f)dfids 
Jo Jm 

^/T{f)d^ids < VUv Vt [ ^/T(f)d^i. 

Jm 



We thus obtain the following basic inequality: for / G Cq{M), 

\\Ptf - /iili(m) <V6di^vi wVnDhHM), t > 0. 



;ii.82) 



Suppose now that C M is a bounded Caccioppoli set. But then, G BVn{Q), for any 
bounded open set D E. It is easy to see (see e.g. the proof of Lemma 2.5 in [2Tj) that 
MarnilE]^) = Var7^(l£;;M), and therefore 1e S BVh(M). By Proposition 111.21 there exists a 
sequence {fn}neN in C^(M) satisfying (i) and (ii). Applying (jll.82p to /„ we obtain 



\Ptfn - /nllLi(M) < Vt || ||li 



(M) 



6di^ Vt Var-H (/„ , M) , n G N. 
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Letting n — > oo in this inequality, we conclude 

||Pil£; - 1£;||li(m) < ^^^^ VarnilE, 
Observe now that, using Pfl = 1, we have 



6du Vi Pn{E;M), t > 0. 



On the other hand. 



We thus obtain 



\PtlE - li?||Li(M) = 2 (^fi{E) - j^PtlEdf?j 
[ PtlEdti= [ {Pt/2lEf dfi. 

Je Jm 



\\PtlE-lE\\L^iM) = '^(KE)- [ {Pt/2lE)^dfi 

\ Jm 

We now observe that (j9.64p and the assumption (1) imply 



p{x,x,t)< <-^, xGM,t>0, 



where C4 = C^^C{d, k, p2)- This gives 



{Pt/2lErdfi< I p{x,y,t/2ydfi{y)] dfi{x) 



C 

^p{x,x,t)^dfi{x)'j <-^ii{Ef. 



Combining these equations we reach the conclusion 



KE) < Vt Pn{E;M) + t > 0. 

Now the absolute minimum of the function g{t) = At"' + Bt^^, t > 0, where A, B,a,P > 0, is 
given by 



Applying this observation with a = ^ , /3 = -5- , we conclude 



13 a 



with 



KE)^ < C^Pn{E 
C^ = {l + D)^ 



2D 



The latter inequality proves (3). 
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The proof that (3) is equivalent to (4) fohows the classical ideas of Fleming-Rishel and Maz'ya, 
and it is based on a generalization of Federer's co-area formula for the space BVh, see for 
instance [3T]. 



Finally, we show that (4) (1). We adapt an idea in [BD] (see Theorem 3.1.5 on p. 58). In 
what follows we let = 15/(1) — 1). Let p,q £ (0, oo) and < ^ < 1 be such that 

I _e 1-d 

p 17 q ' 

Holder inequality, combined with assumption (4), gives for any / G Lipd{M) with compact 
support 



Ilp(M) < II/IIl''(M)II/IIL(M) - (C3||\/r(/)||il(M)j ll/llig(M)- 

For any a; G M and r > we now let /(y) = (r — d(y, Clearly such / G Lipfi{M.) and 
supp / = B{x,r). Since with this choice ||\/r(/)llLi(M) — /^(^(^i '''))^) the above inequality 



implies 



which, noting that + = , we can rewrite as follows 

1 \ 



t^{B{x,r))>[^J M^(x,^))VP^ 



2' 



where we have let a = rJ?n ■ Notice that < a < 1. Iterating the latter inequality we find 



D+ep 

^(i?(x,r))> r'^^U-'2-'^^U^^')-' ^,{B{x,^)T\ G N. 



^From Theorem 12.121 for any x G M there exist constants C{x),R{x) > such that with 
Q{x) = log2 C{x) one has 

niBix, tr)) > C(x)-4'5(^V(^(2;, r)), 0<t<l,0<r< R{x). 

This estimate implies that 

liminf fiiBix,^))"" >1, j;GM,r>0. 
fc— »oo 2^^ 

Since on the other hand X^j^i '^^ ~ ^'^'^ Z^j^iO ~ 1)'^"' — conclude that 

li{B{x,r)) > (2~'e(^+f')c^^Y a; G M,r > 0. 

This establishes (1), thus completing the proof. 

n 
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12 A sub-Riemannian Bonnet-Myers theorem 



Let (M, g) be a complete, connected Riemannian manifold of dimension d>2. It is well-known 
that if the Ricci tensor of M satisfies the following bound for all V G TM 

Ric(y,y) > (d- pi>0, (12.83) 

then M is compact, with a finite fundamental group, and diam(M) < This is the celebrated 
Myer's theorem, which strengthens Bonnet's theorem. Like the latter, Myer's theorem is usually 
proved by using Jacobi vector fields (see e.g. Theorem 2.12 in [19]). 

A different approach is based on the curvature-dimension inequality 

r2(/, /) > -{Lff + (n - l)pir(/, /), (12.84) 
n 

which one obtains from (j5.44p assuming (jl2.83p . By using ingenious non linear methods based 
on the study of the partial differential equation 

2d 



c{F-'-f) = -Lf,l<p< 



d- 2' 

the inequality ()12.84p implies (for d > 2) the following Sobolev inequality 



d 



{d - 2)B? 



(I Ifl^dA - I fdfi < [ T{f,f)df,, /gCo^(M). (12.85) 
\Jm J Jm Jm 



where fi is the Riemannian measure and. Using a simple iteration procedure, it is deduced from 
(|12.85p that the diameter of M is finite and bounded by , see [45] . This non linear method 
seems difficult to extend in our framework. However, a weak version of the Myers theorem could 
be proven by Bakry in [3j by using linear methods only. This method based on entropy-energy 
inequalities (a strong form of log-sobolev inequalities) can be extended to our framework. 

In this section we establish the following sub-Riemannian Bonnet-Myer's compactness theorem. 
Here again, we assume that (M, d) is complete. 

Theorem 12.1 If there exist constants pi,P2 > and a constant k > such that for every 
smooth function f : M. ^ M: 

nf, f) > Pi^if, f) + P2r^if, f) (12.86) 
T{f,f)<KT{f,f), (12.87) 
then the metric space (M, d) is compact in the metric topology with a Hausdorff dimension less 
than d ^1 + and we have 



diam M < 2^7rJ^^±^ ( i + ^] d. 

V piP2 V 2p2y 



We shall proceed in several steps. Throughout this section, we assume that (jl2.86p and (jl2.87p 
are satisfied. 
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12.1 Global heat kernel bounds 

Our first result is the following large-time exponential decay for the heat kernel. 

Proposition 12.2 Let < < There exist to > and Ci > such that for every 

f G C^{M), / > 0; 

d 

' " ^ < Cie-''^ xeM,t>to. 



dt 



In Pt fix] 



Proof. In Proposition 16.51 we choose 

b{t) = (e""* - e-°^)^, < t < T, 
with (3 > 2 and a > 0. With such choice a simple computation gives, 

7(t) = - f 2pi - a/3 - a/3- - e""^ ( a{p - 1) + 

Keeping in mind that 6(r) = 6'(r) = 0, and that 6(0) = (l-e"'^'^)^, 6'(0) = -a/3(l - e-'^^)^-\ 
we obtain from (16.530 



a/3(l - e-"'^)'3-i 



2p2 



-r(ln Pt/) - (1 - e-°^)^r^(lnPT/) 



(12.^ 



Now, 



[ h'{t)-i{t)dt = - ^ f 2pi - a/3 - a/3— ^ (1 - e""^)^ 

+ fa/3 - a + a/3-^ e-"^(l - e-"'^)^-\ 

41- i V P2y 



rb'itMtf 
Jo 



dt = - - f 2pi - a/3 - a/3-^ (1 - e""^)^ 
16 V P2j 

^2 (2pi - a/3 - a/3^) (a/3 - a + a/3^) 



e-"^(l _ e-°^ ) 



1-^ 



^2 fa/3 - a + a/3^ 
16 1^1 



-2aT, 



1 - e 



-aT\l3-2 



If we choose 



then 



a 



2/91 j02 
^(/02 + /t) 



2pi — a/3 — a/3 — = 0, a/3 — a + a/3 — = 2pi — a, 

P2 P2 
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and we obtain from (jl2.88p : 

<^r(lnP^/) + (1 - e-^)r^(lnP^/) < ii^^lZL^.-.T^ ^^^i 



/3, 

Noting that 2pi - a = j^^^{{l3 - l)p2 + (3k) > 0, and that /3 > 2 implies a < ([12:89]) 
gives in particular the desired lower bound for ^ lnPt/(a;) with v = a. 

The upper bound is more delicate. We fix < r/ = ^g^^^^, and with 7 = 2(3pip2 we now choose 
in (fTMl) 

2pi/52 — 76"'^'^ 76"'''^ 



a = — — r — = r/ 



Clearly, q > provided that T be sufficiently large. This choice gives 

K 76 AC 76 

2pi — a/3 — aj3 — = , a(3 — a + af3 — = 2pi — a . 

We thus have 



Jo 4 [ /02 /3 - 1 



a ) \ . 

P2 



Noting that e'^''-")'^ = e /5(''2+-) ^ 1, and a — > -^r^^ as T ^ 00, we obtain 

-(2pi-Q ) — > — -{2pi 



P2 13-1 P2 P2 l3-l\ P{p2 + k) 

Since by our choice of 7 we have ^ — -j^^ ^2/3i — > 0, it is clear that we have 



o \P2 13 - I \ (3{p2 + k) y y 

provided that T be large enough. We also have 



6'(t)7(t)^cit = -^e-^-^(l-e-°^r^|^(2pi-a , 

+ - e-"^)2e-2('?-")r _ 2^-^(1 - e-°^)(2pi - a - 2£:!^)e-(r?-")T\ 

P2 P2P-I p2 J 

Using our choice of 7 we see that, if we let T ^ 00, the quantity between curly bracket in the 
right-hand side converges to 



(3-2 V /3(P2 + k) ) " /? - 1 /3(/52 + k) 
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This quantity is strictly positive provided that 



2(3 (/3 - l)p2 + ^ 1 ai3-l)p2 + l3n V ^ ^ 



P - 1 /3(p2 + k) /3 - 2 V (3{P2 + 1^) 

and this latter inequality is true, as one recognizes by applying the inequality 2xy < + . 

at 



From these considerations and from (jl2.88p we conclude the desired upper bound for -St In Ptf{x) 



□ 

Proposition 12.3 Let < u < ^]^^^ . There exist to > and C2 > such that for every 
f £ C^{M), with / > 0, 

< < gC.e--d(x,y)^ X,y€M,t> to. 

Ptfiy) 

Proof. If we combine (|12.89p with the upper bound of Proposition 112.2} we obtain that for x € M 
and t > to, 

r(lnPJ)(x) < 026-^'. 

We infer that the function ^(2;) = C^^'^e'^* In which belongs to C°°(M), is such that 

T{u){x) < 1, X G M. From (f2Tri) we obtain that 

1^(2;) — u{y)\ < d{x, y), x,y £ M. 

This implies the sought for conclusion. 

□ 

If we now fix X G M, and denote by p{x,-,t) the heat kernel with singularity at (x,0), then 
according to Proposition 112.21 we obtain for t >to, 



d \np{x, y, t) 



dt 



< Ci exp {-vt) , (12.90) 



with < u < ^^^^ . This shows that lnp(-,-,t) converges when t — > 00. Let us call Inpoo 
this limit. Moreover, from Proposition 112.3] the limit, \npoo{x,-) is a constant C{x). By the 
symmetry property p{x,y,t) = p{y,x,t), so that C(x) actually does not depend on x . We 
deduce from this that the invariant measure fx is finite. We may then as well suppose that // is 
a probability measure, in which C3<SG Poo 

= 1. We assume this from now on. 
We now can prove a global and explicit upper bound for the heat kernel p{x, y, t). 

Proposition 12.4 For 2;,y e M and t > 0, 

1 

p{x,y,t) < 



1 _ g 3{p2 + n) 
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Proof. We apply (112. 89p with /? = 3 and obtain 

T{lnPtf) + (1 - e-°*)r^(lnPt/) < ^ + 3k ^,,,L^ ^^^.gi) 



P2 + 2p2 /02 + K ^t/ 

/2p2 + 3Ky 



+ 



12p2 VP2 + «^/ 1-e 



-at ■ 



where a = -^r^^r^- We deduce 



9i ~ 6 p2 + «^ 1 - e-"* ■ 
By integrating from f to oo, we obtain 

- lnp(x, y,t)>-^(^l + ^^ ln(l - e-"*). 

This gives the desired conclusion. 

□ 



12.2 Diameter bound 



In this subsection we conclude the proof of Theorem 1 1 2 . 1 1 by showing that the diam M is bounded. 
Since we have assumed that (M, d) be complete, this implies that such metric space is compact. 
The idea is to show that the operator L satisfies an entropy-energy inequality. Such inequalities 
have been extensively studied by Bakry in [3] (see chapters 4 and 5). 



To simplify the computations, in what follows we denote by D the number defined in (16.56P 
and we set 

2/01 P2 



a 



3(/02 + n) 



Proposition 12.5 For f G L^(M) such that f^f'^dfi = 1, we have 



flnfdfi<<i> 



nf)dp , 



where 



$(x) = D 

Proof. ^From Proposition 112.4] for every / G L^(M) we have 

WPtfWoo < 

(1-e- 



aD \ aD 



D IIJ Il2- 

at\^ 



Therefore, from Davies theorem (Theorem 2.2.3 in [23j), for / € L^(M) such that Jy^f'^dp, = 1, 
we obtain 

f f^ In f^dn < 2t [ r(/)d/z-Z)ln(l-e-"*) , t > 0. 
Jm Ju 
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By minimizing over t the right-hand side of the above inequahty, we obtain 

flnfdfi<—xlnl-——)+Dln( —- 

a \2x + aD J \ aD 

where x = J-^ T{f)d^. It is now an easy exercise to recognize that the right-hand side of the 
latter inequahty is the same as 

□ 



With Proposition 112.51 in hands, we can finally complete the proof of Theorem ll2.11 
Proposition 12.6 One has 



diam M < 2V2\ —n = 2V3tt\ / ( 1 + ^'^ 



.. , _ .d. 

a V PiP2 V 2p2, 

Proof. The function <I> that appears in the Proposition 112.51 eniovs the following properties: 

• ^'{x)/x^^'^ and ^{x)/x^^'^ are integrable on (0, oo); 

• $ is concave; 

• h lo"^ ^dx = /+°° ^dx = -2 /+°° V^'^"{x)dx < +00. 



We can therefore apply the beautiful Theorem 5.4 in [3j to deduce that the diameter of M is 
finite and 



r+oo 

diamM<-2 ^/x^"{x)dx. 
Jo 



Since = — x(^2x+aD) ' ^ routine calculation shows 



-2 V^^"ix)dx = 2^/37r (i + ^)d. 

Jo V ^1/^2 V 2p2y 



□ 



Remark 12.7 The constant 2\/37rW £2±i: m _|_ ^ \ ig ^lot sharp. For instance, in the Rie 



pip2 \ 2p2 

mannian case k = p2 = 0, we obtain 



r- / d 
diam M < 2V37r W— , 

V Pi 

whereas it is known from the Bonnet- My er's theorem that 



diam M< vr^ 



Pi 
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12.3 Dimension bound 

We now turn to an upper bound for the Hausdorff dimension of the compact metric space (M, d). 
Proposition 12.8 The Hausdorff dimension of the metric space (M, d) is less than D given by 

Proof. Let us recall that, from our assumptions, 

TrM = W(x)eV(x), xEM, 

where 

= span{Xi(x), , x G M, 

and 

V(x) = span{Zmnix), 1 < m,n < [).} , 

We moreover assumed that 

dimW(x) =d, xGM, 

and this implies that also dimV(2;) = dimM — d, is independent of x G M. From Theorem 2 
in [50] we deduce that the Hausdorff dimension of the compact metric space (M, d) is equal to 
dimHaus(M) = d + 2dimV(x). Moreover, from [llj and [68] (see also Chapter 3 in [8]), there 
exists a smooth and positive function m on M such that 

limt^p(x,x,t) = m(x). 
Prom the bound (|6.57|) we conclude that 

dimHaus(M) < D. 

□ 

12.4 Isoperimetric bounds and Poincare inequality 

We recall our assumption, following Proposition 112.3] that ^(M) = 1. Also, let D be defined by 
(|6.56|) . With this in hands, we can now proceed as in section [TTl 

Proposition 12.9 Let E cM. be a Caccioppoli set. We have 



f,{E){l-f,iE))<lD.^±^Pn{E,M). 
2 V dpip2 

Proof. We proceed exactly as in the proof of Proposition 111.41 to obtain from (jl2.9ip the 
inequalities 



\\Vmf)\\^<D 



PiP2 e 



-at 



d{K + p2) Vl - e-"* 
and 

11/ - PtfWi < 3D,/^±^Vl-e-"*||/||i. 
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Combining this with Proposition 112.41 gives 



We conclude by letting t — > +00. 

The previous isoperimetric inequality leads to the following Poincare inequality. 
Proposition 12.10 Let f G C°°(M), then 

inf / 1/ - c|d^ < QdJ^^ I yf(7)rf^- 
cgkJm V <ipiP2 Jm 

Proof. Let m be a median for /, that is 

m(/ > m-) > ^, < "i) > ^. 

Set 

= max(/ — m, 0), f~ = — min(/ — m, 0) 
so that f — m = — f~ . We have 



/ \f-m\dfi 
Jm 



and thus 

/* /"H-oo /"H-oo 

\f-m\dfi= fi{f+>t)dt+ fi{r>t)dt. 

Jm Jo Jo 

Observe that for every t > 0, 



from Proposition 112.91 we obtain 



>t)< 3D^I^^±^Pn{{f+ > t},M), 



and 



This gives 



Pif- >t)< 3DJ'^±^Pn{{f- > t},M). 

\f-m\dp<6D.[^P^([ 7fCr)d/.+ / VnF)dp 
V dpip2 \Jm Jm 

Observing that y/T{f+) + yfCP) = \/r(/+ + /~), completes the proof. 
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12.5 A Lichnerowicz type theorem 

A well-known theorem of Lichnerowicz asserts that on a d-dimensional complete Riemannian 
manifold whose Ricci curvature is bounded below by a non negative constant p, then the first 
eigenvalue of the Laplace-Beltrami operator is bounded below by ■£^. In this section, we provide 
a similar theorem for our operator L. Let us observe that in [33J, Greenleaf obtained a similar 
result for the sub-Laplacian on a CR manifold. 

Proposition 12.11 The first non zero eigenvalue Ai of —L satisfies the estimate 

P1P2 



Ai > 



Let us now recall that 



Proof. Let / : M ^ M be an eigenfunction corresponding to the eigenvalue — Ai. ^From our 
assumptions, 

r2(/, /) + i^rf (/, /) > i(L/)2 + (pi - ^) r(/, /) + p^r^if, /). 

By integrating this inequality on the manifold M, we obtain 

/ V2{fJ)dp + v f TUf,f)dfi>] [ {Lffdp+(pr--) [ T{fJ)dp + p2 I T^{f,f)dp. 
Jm Jm " jm ^ Jm Jm 

T{f,f) = ^{Lif^)-2fLf) 
T2{f,f) = ^[LTif,f)-2T{f,Lf)], 

rf (/, f) = l [LT^if, f) - 2T^{f, Lf)] . 
Therefore, by using Lf = —Xif and integrating by parts in the above inequality, we find 

f A? - 4 + ^ - piAi) / f'dp>{p2-.Xi) [ T^ifJ)dp. 
\ a 1^ J Jm Jm 

By choosing = we obtain the desired inequality 

□ 

Remark Vl.Vl We note that when k = 0, we recover the classical theorem of Lichnerowicz. 



and 
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